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ABSTRACT 


Itie  differential  equations  and  boundary  conditions  describing  the  behavior 
of  a finitely  deformable,  heat-conducting  composite  material  are  derived  by  moans 
of  a systematic  application  of  the  laws  of  continuum  mechanics  to  a well-defined 
macroscopic  model  consisting  of  interpenetrating  solid  continua.  Each  continuum 
represents  one  identifiable  constituent  of  the  N-constituent  composite.  The 
influence  of  viscous  dissipation  is  included  in  the  general  treatment.  Altliough 
the  motion  of  the  combined  composite  continuum  may  be  arbitrarily  large,  the 
relative  displacement  of  the  individual  constituents  is  required  to  be  infinites- 


in  the  absence  of  heat  conduction  and  viscosity  is  exhibited  in  detail  for  the 
case  of  the  two-constituent  composite,  Ihe  linear  equations  are  written  for 
both  the  isctroi  '.c  and  trcuisversely  isotropic  material  symmetries.  Plane  wave 
solutions  in  tjie  isotropic  case  reveal  the  existence  of  high-frequency  (optical 
type)  branches  as  well  as  the  ordinary  low-frequency  (acoustic  type)  branches, 
and  all  waves  are  dispersive.  For  the  linear  isotropic  equations  both  static 
and  dynamic  potential  representations  are  obtained,  each  of  which  is  shown  to 
be  complete.  The  solutions  for  both  the  concentrated  ordinary  body  force  and 
relative  body  force  are  obtained  from  the  static  potential  representation. 
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1 . Introduction 

A composite  material  is  composed  of  a i umber  of  distinct  identifiable  con- 
stituents with  different  physical  properties,  which  are  combined  to  form  a single 
solid.  In  particular,  fiber-reinforced  composites  consist  of  fiber  reinforcement 
imbedded  in  a matrix.  In  order  to  obtain  certain  types  of  information,  c.g.,  the 

actual  bond  stress  between  the  reinforcement  and  matrix,  the  individual  elements 

. 1 2 

of  the  composite  must  be  treated  as  separate  entities  * . Nevertheless,  it  must 

be  remembered  that  it  would  be  virtually  impossible  to  consider  an  external 

boundary  of  a composite  while  considering  the  individual  elements  in  completely 

separate  detail.  On  the  other  hand,  in  order  to  obtain  certain  other  types  of 

information,  e.g.,  the  effective  elastic  constants  of  the  single  composite  ma- 

3 4 

terial,  some  sort  of  single  continuum  model  can  readily  be  employed  ’ . However, 
there  are  numerous  other  situations,  e.g.,  the  wave  velocity  dispersion  induced 
by  the  fiber  reinforcement  in  the  matrix  and  possible  resonances  at  whicli  the 
reinforcement  might  separate  from  the  matrix,  for  which  neither  of  the  aforemen- 
tioned approaches  can  adequately  or  conveniently  account  but  a different  model 
somewhere  between  the  two  can  account.  Since  composite  materials  composed  of 
reinforcement  spaced  uniformly  densely  in  the  matrix  can  be  modeled  as  interpene- 
trating solid  continua  and  such  a model  can  conveniently  account  tor  mucli  of  the 
physical  phenomena  that  neither  of  the  aforementioned  approaches  can,  we  employ 
this  model  in  obtaining  a description  of  composite  materials.  At  tnis  point  v/c 
note  that  for  the  intej.penetrat.ing  solid  continuum  model  to  be  valid  a charac- 
teristic length,  such  as  a wavelength,  must  be  large  compared  with,  say,  the 
spacing  of  the  fiber  reinforcement,  in  much  the  same  way  that  the  wavelength  of 
an  elastic  wave  must  be  large  compared  with  a lattice  spacing  for  the  ordinary 


elastic  continuum  description  to  be  valid. 


2, 


The  interpenetrating  solid  continuum  model,  vhich  is  clorjcly  related  to  the 
model  of  fluid  mixtures^,  has  been  employed  in  the  descripi  ion  of  a variety  of 
physical  phenomena  such  as,  e.g.,  certain  types  of  magnetoelastic  interaction^, 
electroelastic  interaction^  and  the  interaction  of  the  electromagnetic  field 
with  deformaVjle  insulators  . In  this  latter  case  in  order  to  consider  ionic 
polarization  resonances,  the  model  consisted  of  two  interpenetrating  continua  in 
which  the  motion  of  the  center  of  mass  of  the  two  continua  was  finite  but  the 
relative  motion  of  each  of  the  continua  with  respect  to  the  center  of  mass  was 
infinitesimal.  It  was  felt  that  in  order  for  the  description  to  be  physically 
meaningful,  the  relative  displacement  of  the  two  continua  had  to  bo  infinitesimal, 
or  else  the  solid  would  rupture.  The  application  of  this  model  to  the  description 
of  material  composites  should  be  obvious.  Indeed,  the  idea  of  employing  inter- 
penetrating continua  as  a model  of  composite  materials  has  already  been  intro- 
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duced  by  Bedford  and  Stern  > . However,  there  are  a nvenber  of  fundamental  dif- 

ferences between  the  approach  of  Bedford  and  Stern  and  that  employed  in  Ref.  8. 

For  one  thing,  Bedford  and  Stern  assu.tie  independent  finite  motions  of  each  con- 
stituent while  in  Ref, a,  as  already  stated,  although  the  motion  of  the  center 
of  mass  is  taken  to  be  finite,  the  relative  motion  of  the  separate  constituents 
is  taken  to  be  infinitesimal.  It  i.s  felt  t’nat  Liie  jpi'occwtiirs  ^TTijvlcycd  in  Ref  ,9 
is  physically  unrealistic  because  the  solid  composite  would  rupture  long  before 
the  relative  displacements  became  large.  Secondly,  in  Ref. 9 a conservation  of 
energy  relation  is  written  separately  for  each  constituent  while  in  Ref.  8 a 
single  conservation  of  energy  relation  is  written  for  the  entire  composite. 

When  a separate  conservation  of  energy  relation  is  written  for  each  constituent, 
the  energy  of  interaction  between  the  constituents  is  not  included  in  the  defi- 
nition of  the  stored  energy  density  nor  do  the  associated  rate  terms  appear  in 
che  formal  expression  for  the  first  law  of  thermodynamics.  As  a consequence. 


3. 


in  Ref. 9 different  temperatures  and  entropy  densities  are  defined  and  a separate 
rate  of  entropy  production  inequality  is  postulated  for  each  constituent,  whereas 
in  Ref, 8 one  energy  density,  one  temperature  and  one  entropy  density  are  defined 
and  one  rate  of  entropy  production  inequality  is  employed  in  the  usual  manner. 

As  a result,  in  Ref.  9 with  the  exception  of  the  defined  volumetric  interaction 
terms  which  are  taken  to  depend  on  tlio  constitutive  variables  associated  with 
all  the  constituents,  tiie  oUiei.  dependent  constitutive  variables  for  each  con- 
stituent are  taken  to  depend  on  the  constitutive  variables  associated  with  that 
constituent  only.  On  the  other  hand,  in  Ref, 8 the  resulting  single  thermodynamic 
equation  for  the  combined  continuum  takes  a form  that  indicates  that  all  dependent 
constitutive  variables,  including  the  relative  stresses  associated  wi.tli  all  the 
different  combinations  of  constituents,  should  depend  on  the  constitutive 

variables  associated  with  all  the  combinations  or  constituents,  clearly,  when 

10 

the  theory  of  Bedford  and  Stern  is  linearized  , the  aforementioned  physical  ob- 
jection concerning  tlie  large  relative  motion  of  the  constituents  vanishes. 

However,  it  should  be  clear  from  the  above  discussion  that  certain  intrinsic 
differences  in  the  descriptions  remain. 

In  this  paper  the  above  discussed  model  of  interpenetrating  solid  continue 
is  applied  in  obtaining  a description  of  a tliree-constituent  composite 
material.  In  the  two-constituent  case  the  model  is  identical  to  the  one  em- 
ployed in  Ref.  8 provided  the  electronic  charge  and  spin  continue  are  omitted 
and  the  ionic  charge  is  ignored.  In  the  three-constituent  case  the  model  is  a 
straightforward  generalization  of  the  two-constituent  case,  and  from  tliere  the 
generalization  to  N constituents  is  obvious.  For  obvious  reasons,  tine  general 
equations  are  determined  only  in  the  three-and  N-constituent  cases.  The  pro- 
cedure employed  in  obtainiriy  the  description  is  exactly  the  same  as  in  Ref.  8 , 
but  in  the  absence  of  the  electromagnetic  field.  However,  in  this  treatment 
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simple  Kelvin- type  viscous  dissipation  is  included.  As  already  indicated,  tlie 
motion  o£  tlic  center  of  mass  of  any  point  of  the  entire  composite  may  be  finite 
while  tlic  relative  displacement  of  any  constituent  from  tlie  mass  center  must  be 
infinitesimal.  Each  constituent  interacts  witli  nciglilx>ring  elements  of  the  same 
constituent  by  moans  of  a traction  vector  associated  with  tliat  constituent.  In 
addition,  each  constituent  interacts  with  all  other  constituents  at  that  point 
by  means  of  volumetric  interaction  forces  and  couples,  both  of  which  are  equal 
and  opposite  in  pairs. 

Tlie  application  of  the  appropriate  equations  of  balance  of  mass  and  momentum 
to  the  respective  continua  yields  the  material  equations  of  motion,  which,  as 
usual,  constitutes  an  underdetermined  system.  Tlie  application  of  tlic  eqviation 
of  tlie  conservation  of  energy  to  the  combined  material  continuum  results  in  the 
first  law  of  thermodynamics,  which,  with  the  aid  of  tiuj  seuonu  law  of  thermo- 
dynamics^^^ and  the  principle  of  material  objectivity^^’ , enables  the  deter- 
mination of  the  constitutive  equations  of  our  nonlinear  description  of  cemposiuo 
materials.  These  constitutive  equations  along  with  tlic  aforementioned  equations 
of  motion  and  tlie  thermodynamic  dissipation  equation  result  in  a properly  de- 
termined system,  which  can  readily  be  reduced  to  (3N-t-  1)  equations  in  (3N+  1) 
dependent  variables.  In  order  to  complete  the  system  of  equations,  jump  (or 
boundary)  conditions  across  moving,  not  necessarily  material,  surfaces  of  dis- 
continuity are  determined  from  the  appropriate  integral  forms  of  the  field 
equations.  It  should  be  mentioned  at  this  point  that  the  resulting  system  of 
nonlinear  differential  equations  and  boundary  conditions  is  considered  to  be 
valid  for  fiber  reinforced  type  composites,  as  well  as  some  other  types,  but 
not  for  laminated  composites  unless  the  thickness  of  a lamina  is  small  compared 
to  a critical  dimension  or  wavelang’-h.  It  is  felt  Jiat  the  equations  provide 


a reasonable  description  of  such  materials  as,  say,  fiber  reinforced  rubber. 


If  in  a iihtir  reinforced  coint'osite  Uic  fiber  in  not  continuous  (chopped  fiber), 
the  description  is  simplified  by  nrglectiny  tlie  traction  in  the  constituent  of 
the  model  representing  tlx  discontinuous  (chopped)  fiber  reinforcement  and  in- 
cluding only  the  volumetric  interaction  between  this  constituent  and  the  matrix. 

Since  tlio  resulting  general  noiilinear  equations  are  relatively  intractable 
fer  tl'.e  treatment  of  many  problems  in  their  natural  form,  the  linear  version  of 
the  equations  is  extracted  from  tlie  general  one.  Tlrcse  linear  equations  are 
specialized  to  the  important  case  of  a transversely  isotropic  material,  which 
occurs  very  frequently  in  continuous  fiber  reinforced  composites,  and  the  fully 
isotropic  case.  The  static  potentials  analogous  to  the  Boussinesq  Papkovitch 
potentials  of  c.iassical  elasticity  are  obtained  from  tlie  linear  equations  for 
the  two  constituent  isotropic  composite.  From  this  static  potential  roprcscnl  a- 
tion  tlie  solutions  for  concentrated  forces  in  the  infinite  two  constituent  iso- 
tropic composite  are  obtained.  The  dynamic  potentials  analogous  to  the  Lame 
potentials  of  classical  elasticity  are  obtained  from  the  linear  equations  for 
the  two  constituent  isotropic  composite.  Comi;loteness  is  established  in  both 
cases.  Plane  wave  solutions  of  the  linear  equations  for  the  two  constituent 
composite  are  presented  for  the  isotropic  case.  The  solutions  reveal  the  exi-t- 
ence  of  tv.'c  sets  of  waves,  higher  (ir.  frequency)  ones  and  lower  ones,  as  expected. 
The  lower  ones  are  dispersiv-e  and  approach  the  non-dispersive  velocity  of 
classical  linear  elasticity  from  below  as  the  wavotiumber  anov-iarhcc  cc.o  and  the 
upiper  ones,  which  are  highly  dispersive,  have  non-zero  cul.-oft  frequencies  cor- 
responding to  the  aforementioned  resonances  at  which  the  reinforcement  might 
separate  from  the  matrix. 

Since  the  defined  material  constants  are  not  known  for  any  two  constituent 
composite,  it  is  suggested  that  plane  wave  measurements  be  made  and  correlated 
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'with  Uie  above-mentioned  plane  wave  solutioiio  in  order  to  obtain  the  defined 

linear  elastic  constants  of  tlie  two  constituent  composite,  in  much  tlic  same 

15-17 

manner  as  in  anisotropic  elastic  (or  piezoelectric)  materials  . Hie  theory 

can  then  readily  be  checked  by  comparing  calculations  witli  measurements  in  addi- 
tional redundant  directions.  Moreover,  when  known  material  constants  arc 
available,  such  things  as  surface  wave  velocity  dispersion  can  be  calculated  and 
compared  with  ir.easurement.  In  addition,  systematic  dispersion  information  can 
be  used  for  the  non-destructive  testing  of  fiber  reinforced  composite  materials. 

An  analysis  of  surface  waves  in  a two  constituent  composite  has  been  made  and 
some  calculations  have  been  performed  when  the  fiber  reinforced  composite  material 
is  simplified  sufficiently  that  some  of  the  constants  can  be  estimated  from  t)ic 
knev.-p  elastic  constants  of  the  individual  constituents  of  the  composite.  This 
investigation  will  be  reported  in  a forthcoming  work, 

Xn  closing  tlie  Introduction  we  note  tJiat  more  general  single  continuum 

18“  2 1 

theories,  commonly  referred  to  as  microstructure  tlieories  , can  be  and, 

22 

indeed,  have  been  applied  to  certain  composites  to  describe  some  of  tlic  phe- 
nomena that  the  model  of  interpenetrating  solid  continua  desr->-ibo3.  However, 
the  resulting  equations  are  quite  different  from  tliose  presented  here  and  wc 
find  it  difficult  to  idencify  the  physical  meaning  of  Uie  microstructure 
variables  with  any  degree  of  certainty. 

2 , The  Interacting  Continua 

As  indicated  in  the  Introduction,  the  macroscopic  model  we  first  consider 
consists  essentially  of  three  distinct  interpenetrating  solid  continua.  initi- 
ally, all  continua  occupy  the  same  regiori  of  space  and,  hence,  liavc  the  same 
material  coordinates  X . The  motion  of  the  center  of  mass  of  tlic  combined 

Li 
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continuum  is  described  by  the  mapping 


■ 


(2.1) 


which  is  one-to-c;ie  and  differentiable  as  often  ao  required.  In  (2.1)  the 
denote  the  spatial  (or  present)  coordinates  aixi  X , the  material  (or  reference) 
coordinates  of  tlie  center  of  mass  ar«d  t denotes  the  time.  We  consistently  use 
the  convention  that  capital  indices  denote  the  Cartesian  components  of  X and 
lower  case  indices,  the  Cartesiem  components  of  y.  Thus,  X and  y denote  th<‘ 
initial  position  of  all  material  points  and  the  center  of  mass  of  thic  combined 
continuum,  respectively.  Both  dyadic  and  - trtcsian  tensor  notation  are  used 
interchangeably.  A comma  followed  hy  an  index  denotes  partial  differentiation 
;.'ith  respect  to  a coordinate,  i.e.. 


yi,L  " *K,j  “ 

and  the  summation  convention  for  repeated  tensor  inoiccs  is  employed . The 
superscripts  1,2,3  are  used  to  denote  the  respective  continua.  Since  each 
continuum  possess«=‘S  a positive  reference  mass  density  r (n  “ 1,  2,  3 ) and  initially 
occupies  the  seune  region  of  space,  we  have 

P - = c.  + p_  + r_  , (2.3) 

W U V 


whe-e  is  the  total  reference  mass  density  of  the  combined  continuum. 

In  a (finite)  motion  each  continuum  is  permitted  to  displace  with  respect 
to  the  center  of  mass  of  the  combined  continuum  by  infinitesimal  displacement 
fields  • A schematic  diag'-am  indicating  the  motion  of  tlie  model 

appears  in  Fig.l.  The  infinitesimal  displacement  fields  are  regarded  as 

functions  of  y and  t.  Since  tlie  w^''^  are  infinitesimal  and 

<v 

y = y (y,t)  , 


(2.4) 


8. 


and  the  determinant  of  a matrix  product  is  equal  to  the  product  of  the  deter- 
minants, we  have 

=»  V(1  + V • Rs  V , (2.5) 

where  is  the  present  volume  of  the  nth  constituent,  V is  the  present  volume 

of  the  center  of  mass  and 


V « JV  , 
o ’ 


(2.6) 


where  V is  the  reference  volume  of  the  material  and  as  usual 
o 

J = det  y . 


(2.7) 


Inasmuch  as  mass  is  conserved  separately  for  each  constituent,  from  (2.5)  and 
(2.6)  we  have 

/"’j-oi"',  (2.8) 


which  enables  us  to  write 


(1)  ^ „(2)  ^ (3) 

P * P + P P ^ 


(2.9) 


where 


PJ  = Po  . 


(2  ’0) 


and  p is  the  total  present  mass  density  of  the  combined  continuum.  Since  y has 
been  defined  as  the  center  of  mass  of  the  combined  continuum,  we  may  write 


,<1) 


(y  + w^^^)p^^^dV  + r (y  + w^^^)p^'^  dV  + 
" :,(2)  ~ 


i (y  t w(^^p<3-'dV=  f tp^^^ydV, 

i(3)  - - i 


(2.11) 


and  by  virtue  of  (2,5),  we  have 


+ p">w">  . 0. 


9. 


In  addition,  because  mass  is  conserved  separately  for  each  continuum,  wo 
further  obtain 


,<i)  !L 

^ dt 


(1) 


(2)  ~ 


(3) 


<3) 


+ p ■ ■ 4 p-  ' = 0 , (2.13) 


dt 


where  d/dt  denotes  the  material  time  derivative . 

The  interpenetrating  continua  interact  with  each  other  by  means  of  defined 

local  equal  and  opposite  force  fields  F =-F  , F =-F  . F =-F  , 

which  are  located  at  tie  position  y,  wliere  the  first  superscript  denotes  the 

continuum  being  acted  on  and  the  second,  tl\e  continuum  prcxJucing  the  action,  and 

defined  equal  and  opposite  local  r.aterial  couples  C =-C  ,C  =-C 
L 23  L 32 

C " = - C . Each  continuum  interacts  with  neighboring  elements  of  the  sam<^ 
continuum  by  mear;S  of  a traction  force  per  unit  area  t^^^,  t^^\  acting 

across  the  surface  of  separation.  Schematic  diagrams  illustrating  the  above- 
mentioned  interaction  in  the  model  are  shown  in  Figs. 2-4. 


3 . The  Equations  of  Balance 

In  view  of  the  discussion  in  Section  2,  the  equations  of  the  conservation 
of  mass  for  the  different  continua  may  be  written  in  the  form 


dt 


(1) 


= n 


A.  ^ 

dt  j 
V 


dV  ^ 0 


_d_  : 

d t J 
V 


.n 


(3) 


dV  = 0 , 


(3.1) 


Where  V is  an  arbitrary  element  of  volume  for  which  each  of  the  cdnt.i'ma  has 
the  same  material  coordinates.  From  (2.3)  and  (3.]),  wi}  obiain  tlic  cquathm 
of  the  conservation  of  mass  for  the  combined  continuum  in  the  form 


A r 

dt  J 
V 


p dV  = 0 . 


(3,2) 


The  equations  of  the  conservation  of  linear  momentum  for  each  oi  Uic  ttire(i 
continua  are,  respectively 
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dw'  \ 

j t<»as.  J (V^V^av-A|  -i^)av  , 


(3.3) 


<*  dw'  \ 

j t^^^ds+  J j (V^-^ Tt  J 


r\  N*w  V 

Jt^^^ds+  J p^^^f‘'^dv+  J (V^  + ^ I ■" 

V 


(3.5) 


wherft  V dy/dt.  The  equations  of  the  conservation  of  angular  momentum  for  each 
of  the  three  contir.ua  axe,  respectively 

r „ . „<»,xt<^>as  . j (y  . w<»)xp<^>£">av.  J <V^*V^av 


. I ..X,  av  ^ ..1  r („  4.  ^ y n v + 

V V 


dt  J 
V 


\~  dt 


(3.6) 


f (y  + w<^^Xt^2^dS+  J (;^  + J + 

V <v  ^ 


dw 

+ [ yx  (V^+ V^)  dV=  r (y X V + -^^;  dv  , (3.7) 


J A.  I \ 

O '\f  <V 

V 


- /-3\  i-i\  r r t 31  L-32.  ... 

1 (y  + w'"')Xt'-"dS+  I (y  + W'-')XP’  'f'  ■ dVS-  J ( C + J.  ; av 

s V V 


- I (V^  V6  av  = IJ  <y  X p ">(v  . :=_; 


.4  (3) 
dw 


(3 .8) 


Application  of  (3.3)  - (3.5)  to  an  elementary  tetrahedron  in  the  usual 
manner  yields  the  definition  of  the  respective  stress  tensors  of  each  cf  the 


three  continua,  thus 


,(!).„.  ,(1),  .„  • T«>,  t»>  =„  • T»>. 

> i.  ^ .X.  ^ -V/  ' ~ ~ -X, 


(3.9) 


11, 


Substituting  from  (3.9)  into  (3.3)-  (3.5),  respectively,  taking  the  material 
time  derivatives,  using  (3.1),  the  divergence  theorem  and  the  arbitrariness  of  V, 

we  obtain  2 (1) 

<3  w 


V • T 

tst 


dv 

/N# 

.2  (1) 
d w 

- P 

dt 

P 

, 2 

(2) 

dv 

<3t  , 

d2w<2^ 

- P 

dt  “ 

P 

dt^ 

s . T<^>.  57-  p"’ 

r\t  Q ll 


d2„<3> 


bp2J.^  Lp23  ^ 0 ^ ;3.11) 


, V^+V^  = 0,  (3.12) 


which  are  tlie  stress  equatiorj  of  motion  of  tlie  three  continue,  and  where 
V = e.  d/Sy.  and  e.  is  a unit  base  vector  in  the  ith  Cartesian  direction. 

#v>  X 'vi 

Substituting  from  (3.9)  into  (3.6)  - (3.8),  respectively,  taking  the  material 
time  derivatives,  using  (3.1),  the  divergence  theorem,  (3.10)-  (3.12)  and  the 
arbitrariness  of  V,  we  obtain 

(1)^  ,(1)^(1), 

®X®£ij  ij  ^®r£kj'  k ^ij  ^i  ~ ^ ~ 

..(1)  . .U)  1 „ s.  .(1)  0 , ^3.1: 


- w'  X p 


(3.13) 


(2)  (2)  (2).  (2)  (2)^(2)  . L 21  . L 23 

e .e , T + e .e  . (w  T..).-tw  Xp  f +C  + C 
®£  Xij  i:  ~£  Xkg  k 13  ,x  -V  ^ ~ ~ 

^ 2 (2) 

/ox  /OX  ° y. 

- w'"'  X p'"'  - W'*"xp""  ; ^0  , 

at. 

(3)  , (3)  (3),  4.„(3)  V 4. 

~£  £13  13  £^3  k 13  , 1 " ~ 

-.»>xp">g-w'^'xp<=”^.0, 

~ dt  - ^^2 


(3.14) 


(3.15) 


which  constitute  the  equations  of  the  conservation  of  exngular  momentum  of  each 
of  the  three  respective  continua. 

Adding  (3.10)-  (3.12),  we  obtain 


7 • T + pf  = p —■  , 

^ «%»  at 


(3.16) 
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which  are  the  stress  equations  of  motion  of  the  combined  continuum,  and 

T - T + T + T (3. 17) 

(3.18) 

where  t is  the  total  mechanical  stress  tensor  emd  f is  the  total  body  force  per 


(l).(l)  ^ «(2)  (2)  ^ „(3),(3) 
pf“p  f +p  f +p  f 


unit  mass.  Now  let  us  define  the  constants  r^^^  and  r^^^  by 


(1)  (l)/„(3)  ^(2)  ^(2).  (3) 

r B p /p  , r » p /p  , 

and  then  the  subtraction  of  r^^^  times  (3,12)  from  (3.10)  yields 

/V  ^ ~ Z 'W  '' 

oti 


vhere 


^(1)  ^ ,(1)  _ ^(l),(3)  _ 

ij  ij  13  - 

V^, 

^ »V  ' 


(3.19) 


(3.20) 

(3.21) 

(3.22) 


-d)  ^ ,(1)  _ ,(3) 

f\d 


(3.23) 


(3.24) 


and  we  have  employed  (2.12),  similarly,  subtracting  r‘  ' times  (3.12)  from  (3.11), 
we  obtain  „ .„. 


.2  (2) 

V . + p<2>f<2>-p<2>  -V 

/V,  i-v»  #V  2 

dt 


f .r^2>  . 0 


where 


(3.25) 


i')  [2)  (2)  (3)  ,(2)  L21  L23  (2)  b 31  L32, 

D.  =T,.'-r  T..  , J =s  F + F -r  (F  + F ). 

ij  3.3  13  ' ~ ~ 

-(2)  ^f(2)_^(3)  ^(2>  .,(2). , (2)j^(2)^  ^(1)^(1)  ^ 

•~\j  /N.-  ' rv»  #%/  'N#  rv  *'.< 


Equations  (3.20)  and  (3.25)  are  called  the  difference  or  relative  equations  of 

motion,  and  ofV,  and  are  tJie  difference  stresses  and  difference 

» 13  ’ 13  -J  ' -3 


13. 


dispJ.aceraents,  respectively.  Adding  (3.13)-  (3.15)  and  employing  (2.12)  (3.17), 

(3.1''}-  (3.21)  and  (3.23)-  (3.26),  ve  obtain 

— i Ir]  13  £kj  k ij  k 1]  ,1 

- ;5<^>).w<2>X(V.D<2^^'2>)=0,  (3.27. 

^v*  •'Srf  »V  ^ ^ 

which  is  the  equation  of  th^  conservation  of  angular  momentum  for  the  combined 
continuum.  Equation  (3.27)  turns  out  to  be  of  considerable  value  and  interest 
when  viscous  type  dissipation  is  considered.  However,  in  the  absence  of  viscous 
dissipation  Eg,  (3.27)  is  a direct  consequence  of  tlie  invariance  of  the  stored 
energy  function  in  a rigid  rotation. 

Although  we  cannot  explicitly  evaluate  each  of  the  defined  couples  of 
interaction  between  the  respective  continue  in  the  description  presented 

here,  we  can  readily  evaluate  the  total  internal  couple  acting  on  each  continuum, 
i.e.,  tlie  where 

Lcdl  , V2._I.p3  l,p2)  ^ I-pi,.  V3  Lpi  . I-pi.  V2  13 

<V  ^ ^ O/  r\/  /%■/  ' •%/  i“v»  ^ ’ 


and  that  is  all  that  is  required  in  this  type  of  description.  Specifically, 

the  may  be  determined  a posteriori  from  (3.17),  (3.21),  (3.26)^  and 

(3.13)  - (3.15).  Similar  statements  hold  in  the  case  of  rhe  defined  forces  of 
L mn 

interaction  F between  the  respective  contjnua,  and  we  can  readily  evaluate 
a posteriori  the  total  internal  force  acting  on  each  continuum  from 

(3.17),  (3.21),  (3. 26  )j^  and  (3 .10)  - (3 ,12) , where 

LfCi)  ^ v^V^  - V^V^  = V^V.  (3.20) 
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4 . Thermodynamic  Considerations 

The  conservation  of  energy  for  the  combined  material  continuum  can  be 
written  in  the  form 


^ j ^ P" 

V S 

,(3) 


t<^>  . (v. 

j dt  j 'v 


dw 


(2) 


\ 

dt 


J + ^r")  - 5 • s]  ^ ! [p • (x  + ir")  * 


(2) 


13) 


■ (x  + %r-)  + "'’’j'"’  • (x  + %-)] 


'where  T is  tlie  kinetic  energy  per  unit  volume^  e is  the  internal  stored  energy 

.(n) 

ri  W 

(n) 


dw 


per  unit  mass,  t 


V + 


dt 


n = 1,2,3,  denote  the  rates  of  working  per 


Unit  area  of  the  mechanical  surface  tractions  acting  in  tiie  three  continua, 
respectively,  n * q is  the  rate  of  efflux  of  heat  per  unit  area  and  • 

, , ~ -V  ~ 

f n I 

dw 


( 


\ 


V 

V + “TT — ) denote  the  rates  of  working  i-er  unit  volume  of  body  forces  acting  in 
at  / 

three  continua,  respectively.  In  order  to  obtein  expressions  for  T,  we  must 
return  to  our  model  of  the  combined  material  continuum. 

From  the  model  of  the  continuum  it  is  clear  tlnat  the  kinetic  energy  per 


unit  volvune  is  of  the  form 


<^\'v 

dw'-^' 

dw'^' 

'V* 

-'l  r n 

+ 

dt 

; U 

dt 

J + P 

(3) 

dw<^\ 

+ p 

dt  J 

“dt  , 

dw 


(2) 


V + 


dt 


-)•( 


dw 


(2) 


v + 


dt  / 


(4,2) 


Expanding  terras  in  (4,2)  and  employing  (2.9),  (2,13)  and  (3.19),  we  obtain 


= i [p  V • v + (1+  r^^^ 

2 L 


dw^^' 


dw^^^  dw 


dt 

(2) 


dt 


dt 


dt 


(1)  (2) 


dt 


dt 


■]■ 


(4.3) 
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Substituting  from  (2,13),  (3.18),  (3,23),  (3.26)  and  (4.3)  into  (4.1),  we  obtain 


d r ri  i <!•)  ^ 

^ ■rpv'v  + ^-p  (1+r  ) 

ut  v'  L2  -N/  ^ 2. 

V 


dt 


dt 


dw<2)  dw<2) 


dw(-) 


dw^^^  dw 


2 ^ dt 

r / (1)  ' ' - (2) 

! (t  * V + d'  ^ ' -T~  + a'  ' • 

\r^  ru  <w  Cl  U 

S 

.(i) 


^ t . pe]  dV  = 


P 

(2) 


dt 


r / n^~ni  (2)~(2) 

+ J vpf  • p i 

V 


n * q ) dS 

(2) 


dt  . 


•)  dV  , 


(4.4) 


where 

t . t'U  + * t<=>  d<»  . t“>  - , 

(4.5) 

and  from  (3.9),  (3.17),  (3.21)  and  (3.26)^,  we  have 


t * n 


T , d 


= n • , d'^^  * n • D 


(2) 


/V  ^ 


(4.6) 
(n) 


and  t is  the  mechanical  traction  vector  of  the  combined  continuum,  and  the  d 
n = 1,  2,  may  be  thought  of  as  the  difference  traction  vectors.  Taking  the 
material  time  derivative  in  (4.4)  and  using  (4.6),  (3.16),  (3.20),  (3.25)  and  employing 


zr\e  aivergence  ■cneuiem^  unc  luaueLxai  i^xmc  j.vaua.v»=io  wj. 
the  arbitrariness  of  V,  we  obtain 

de 


=s  T . .V 


dt  ij  j 


, , d»'U^ 

..  3 ^ 4.  T^'  ^ > 

,i  ^ij  ^ dt  /,i  ij  V~<TE 


- . 


dw 


(1) 


:(2) 


dw 


(2) 


dt 


dt 


- q, 


(4.7) 


which  is  the  first  law  of  thermodynamics  for  our  combined  continuum. 

We  may  now  introduce  dissipation  by  assuming  that  the  symmetric  part  of 
the  total  stress  tensor  the  two  difference  stress  tensors  and 
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and  the  two  difference  internal  forces  3^  ^ and  3^  ^ may  be  written  as  a sum  of  a 

dissipative  and  a nondissipative  part.  This  is  a restrictive  assumption,  but 

24 

it  is  believed  to  bo  adequate  for  our  purposes  , Thus,  we  write 

S R S ^ D S ^(1)  II  (1)  ^ D (1)  (2)  _ R^{2)  D (2) 

T=*T  + -T,D  = D + D ,D  = d + d , 

' 'Vi  'V' 

^ ^5(^5  , 3<^>  <4.8) 

-VI  rv  ^ ^ *V.  'V-  ' 

and  in  each  case  the  superscript  R indicates  the  nondissipative  (stored  energy) 
portion  and  the  superscript  D,  the  dissipative  portion.  Substituting  from  (4.8) 
into  (3.27)  and  obtaining  the  tensor  form  from  the  vector  form,  wo  obtain 


A R A DA 
T . . = T . . + 1 . , 

13  IJ  13 


(4.9) 


where 


1]  2 V ki  g,k  kg  i,k  13  3 ^ 

ki  g,k  kg  t,k  1 g g 1 J ' 


13 


2 L ki  j,k  kj  i,k  X g g 1 


+ °D<2)^(2)  _ V2)^(2)  _ Dg(2)„(2)  ^ D3(2)^(2)- 

ki  g,k  kg  i,k  x g g x J 


(4.10) 


S A 

Since  T = T + T ’ from  (4.8),  and  (4.9),  we  may  write 

i-w  ^ ^ 

R D S DA 
T = T + T + T , 


(4 .11) 


where 


R R S R A 

T = T + T 


(4.12) 


Equation  (4.11)  is  the  form  we  are  interested  in  employing  because  it  enables 
us  to  obtain  all  the  results  of  interest  to  us  most  readily.  Substituting  from 
(4.8)2  ^ and  (4.11)  into  (4.7)  and  employing  (4.10),  we  obtain 


■ X. 
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p ,v,  . . CM  • - ^ \P—)  ■ 

dt  13  j,i  JL3  \ dt  J ,i-  3 dt  13  \ ot  /,i 


- 

2 


, (2) 
aw . 

D s , 

T Q 

[( 

dt 

ij  ij 

kj  L' 

^ dt  /,k  i,k  ijj 

(1)  1 
T .r  - /o 

dw<2) 

L dt 

pdw<2) 

r 3 

* W . w , . 

1 13  J 

+ 

1^  ‘^’kj 

L \ dt  / jk  i i j J 

(2)  1 

L dt 

- w.  uj.  . 

1 13  J 

1 - 

J 

(4.13) 


vh-sre  and  are  the  rate  of  deformation  and  spin  tensors j respectively. 


which  are  defined  by 


d..=i(v.  ,+v.  .),  u)  . = 4(v.  .-V.  .), 

13  2 3,1  1,3  13  2 3,1  1,3 


(4.14) 


For  the  circumstances  we  have  outlined,  tlie  mathematical  expression  of  tlie 
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second  law  of  thermodynamics  may  be  written  in  the  form 


de  R 
P -T-  ■"  T , .A 
^ dt  13 


.V.  . - ( 3 ) . + 3 ) ■ 

j 3,1  13  \ dt /,i  3 dt  13  \ dt /,l 


1^(2)  ""3  _ clTl 


(4.15) 


where  6 is  the  positive  absolute  temperature  and  1]  is  the  entropy  per  unit  mass. 
From  (4.13)  and  (4.15),  we  have  the  dissipation  equation 

^^d.  . ^ V - -1  - , 

13  ^3  ^3  Lv  dt  /,k  i,k  13J  3 L dt  1 *^13 

du,(2)  .„(2) 

+ V^^r(— i— ) V - - q.  . = pe  “ . (4.16) 

k3  L\  dt  /,k  i,k  13J  3 l dt  1 13J  1,1  dt 


and  the  entropy  inequality  may  be  written  in  the  form 


18. 


(4.17) 


where  T is  the  positive  rate  of  entropy  production.  At  this  point  it  should  he 
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noted  tliat  this  theory  can  readily  be  generalized  * to  account  for  more 

general  functional  constitutive  response  in  Uie  manner  set  forth  in  a previous 
7 

paper  . 


5 , Constitutive  Equations 

Since  we  arc  concerned  with  thermodynamic  processes  for  whicli  both  the 
state  function  equation  (4.15)  and  the  dissipation  equation  (4.16)  are  valid, 
we  may  determine  the  dissipative  constitutive  equations  from  (4.17)  and  the  r.on- 
dissipative  constitutive  equations  from  (4 .15)  which,  by  virtue  of  the  relations 


^M,i  dt 


j,M  ’ \ dt  /,i  Mji  dt  ^ \ dt  /,i  M.i  dt  j,M 


may  be  written  in  the  form 

de  „ d , , FL(1)  , <3  , (1),  K,(2)  d , (2) 

dt  ~ ij  M,i  dt  ^j,M  ij  ^M,i  dt  '^j,M  ^ ij  ^M,i  dt  ,M 


j dt 


dw<2) 

R.>(2)  j 
j dt 


dT| 

p®  d^ 


(5.1) 


Since  the  entropy  inequality  is  of  the  form  shown  in  (4,17),  it  turns  out  to  be 
convenient  to  define  the  thermodynamic  function  \ji  by  the  legendre  transformation 


^ = e - TlQ  . 


(5.2) 
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The  subotitution  of  the  material  time  derivative  of  (5.2)  into  (5.1)  yields 


dw'^’  dw^2) 

Itjd)  ^ j IU(2)  ^ j 

j dt  “ j dt 


•n 

dt  • 


(S.3) 


Since  (5.3)  is  a state  function  equation,  we  must  have 

4 - My  ■ w<2>  . w<2).  e) 


(5.4) 


Substituting  the  material  time  derivative  of  (5.4)  into  (5.3),  we  obtain 

^ d , . . /It.(l)..  S|- 


3,n  ^ j,M^ 


j >M 


St _d^  (1) 

(] 
j 


(%"’  * p ;%t)  ^ -r’  - p(^  * ID  i - ■ 


(5.5) 


3w.' 

3 


since  (5,5)  holds  for  arbitrary  d(y.  „)/dt  d(wf^^)/dt,  d(wf^^/dt,  dwf^Vdt, 

■'i,M  3,M  ' g,^r  ’ 3 


(2) 

dWj  /dt  and  d6/dt,  wo  have 


P,_ 


■ P *♦''3  , V"'  --  H S«/dwj"’  , 

= P , T1  = - p d,/se  . 

M,x  13  g,M 


(5.6) 
(5  7) 
(5.8) 


R R_  ^ 1 ) R_  ( 2 ) 

Solving  (5.6)  , (5.7)  and  (5.8),  for  t , d'  and  D , respectively,  wc  find 


V<“  =py. 

■3,M 


13  "i.Md(y^„)  ' ■'ij  ^ (l)j  ’ “ij  '-"i,M  ^^^(2)j 

" '*  j -M 


9t IS^(2)_„ 

,<!))  ’ 

3,m' 


St 


(5.9) 


Clearly,  \|i  cannot  be  an  arbitrary  function  of  the  variables  shown  in  (5.4) 


13  14 

becau.'e  in  order  to  satisfy  the  principle  of  material  objectivity  ’ 


e and . 
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hence,  i|i  must  be  scalar  invariant  under  rigid  rotations  of  tJie  deformed  body, 
and  any  arbitrary  function  of  the  34  assumed  variables  (11  vectors  and  a scalar 
at  the  point  yj^)  will  not  be  so  invariant.  However,  there  is  a theorem  on 
rotationally  invariant  functions  of  several  vectors  due  to  Cauchy^^,  which  says 
that  iji  may  be  an  arbitrary  single-valued  function  of  the  scalar  products  of  the 
vectors  and  the  determinants  of  their  components  taken  three  at  a time . Appli- 
cation of  this  theorem  shows  that  'ij  is  expressible  as  an  arbitrary  function  of 
66  scalar  products  and  127  determinants  and  0 for  a total  of  194  quantities. 
However,  the  194  quantities  are  not  functionally  independent  and  it  is  relatively 
easy  to  show,  by  using  procedures  similar  to  those  employed  in  Section  6 of 
Ref . 6 that  the  194  variables  are  expressible  in  terms  of  the  31  arguments 
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^k,L  k,M 


V,  w. 
^k.L  k 


(1) 


LM 


. (2) 

yk,L''k,M 


,,<2) 

= ^k,L\ 


(5.10) 


Thug,  we  find  that  i)!  is  invariant  in  a rigid  rotation  if  it  is  a single-valued 
function  of  the  31  arguments  listed  in  (5.10).  Hence  \|;  may  be  reduced  to  the 
form 


* = nr.  0) 


LM  ’ LM 


(5.11) 


in  place  of  tlie  form  shown  in  (5,4),  and  v/herc  we  have  taken  the  liberty  of 

replacing  Green's  deformation  tensor  C , whicli  does  not  vanish  in  tlie  undo- 

KL 


formed  state,  by  the  equivalent  material  strain  tensor  E , which  docs  vanish 

IGj 


in  the  undeformed  state,  and  is  related  to  C by 

KL 
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^KL  ~ 2 ^^KL  ’ 


(5.12) 


From  (5.6),,,  (S.V)^,  (S.S)^  and  (5,9)-  (5.12),  we  obtain 


^..=py.  y.  ^ .py,  ^_w.<">  + py.  -It- 

13  ^ 1,L  ^^a)  3 ^„(2)  3 


aN 


+ py  _^„(l)  + py  -ll-  „(2) 

P^i,L  (1)  j,M  (2) 

LM  LM 


ij  ._(1^  ^ ij 


L 3p(2) 
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9j( 1^(2)  _ 
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r 
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(5.13) 


(5.14) 


(5.15) 


Where  we  have  Introduced  the  conventions  a'l/aE.. ..  = ait/aE,.,  and  it  is  to  be 

JjM  ml 

assumed  that  dE  /dE  = 0 in  differentiating  ijr.  Substituting  from  ( .14)  and 
XL  LK 

(5.15)  into  (5,13)  and  employing  the  chain  rule  of  differentiation , we  obtain 


T, . = PY,  rV 


ij  ^-'i,L-'j,M  aE, 


^ - ^.^2)^(2) 
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V'^wf2)  . 

ki  j,h  ki  3,k 


(5.16) 


Clearly,  the  antisymmetric  part  of  t.  obtained  from  (5.16)  is  identical  witf. 

R ’V 

the  expression  for  given  in  (4.10).  Taus,  even  in  this  rather  conplex 

situation,  tlie  aJitisyTTirnetric  portion  of  the  nondissipativc  part  of  the  stress 
tensor  is  derivable  from  a thermodynamic  state  function  and  has  just  the  value 
required  by  the  conservation  of  angular  momentum. 

Tills  brings  us  to  a consideration  of  the  dissipative  constitutive  equations, 
v/hich  are  obtained  from  tire  entropy  ineq-aality  (4.17)  which  may  be  written  in 
the  form 


^D^(2)  U)  e 

13  13  k3  =k3  3 3 ^3^43  33  ,1 
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where 
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— dwf^Vdt  - 
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3 

3 

1 

,(2) 

'■  i OD-  ■ > 

= dwf^Vdt  - 

w.<2) 

i,k  i: 

3 

3 

1 

~ = dw^(^Vdt  - wf  >o>i,  . <5.18) 

Motivatad  by  (5.17)  we  take  the  dissipative  constitutive  equations  in  the  fonn 

M ft<2)  ft  X 


.(1)  q(1)  A2)  ,,(2) 


D (1)_D  (I)  (1)  a.)  A2)  (2) 

j ~ '^D  > Pk  > ’ K ’ ®,k^  > 

Dr,{2)  _D  (2)  ,,  ^(1)  .(1)  ^(2)  (2) 

^ij  - ^ij  kjL  ’ > ^JL  ’ \ > ^k^ 

D„(2)^D  (2)  (1)  (1)  (2)  (2) 

~ '^j  ^\V  ^kH.  ’ \ ’ \JL  > \ ’ ®,k>  > 

^i  " '^i^'^k£’  ^k^  •’  ^k  ' ^k^  ’ \ ’ ®,k'  •’  i=..l9) 

but  since  the  nondissipative  constitutive  equations  (5,13)-  (5.15)  depend  on 

the  y,  and  0,  there  is  no  .logical  reason  to  exclude 

1 ’ 1 ’ i,M  ' ^ 

33 

them  from  the  dissxpative  constitutive  equations  . Hence,  on  account  of  the 
chain  rule  ot  differentiation,  wc  may  write 
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n (1)  (1)  (2)  (2)  . (1)  (1)  (2)  ^(2)  . 

D^^(2)  (dj^^,  , Pj^  ' ®,M’  ^k,M'  k ' k,M'  k ' k,M’ • 
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(5.20) 


23. 


for  the  general  functional  dependence  of  the  dissipative  constitutive  equations. 

Now,  in  order  that  the  dissipative  portions  of  the  constitutive  equations  satisfy 
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the  principle  of  material  objectivity  ’ , all  variables  in  (5.20)  roust  be 

objective,  i.e.,  they  must  transform  as  tensors  under  time-dependent  proper 

orthogonal  transformations.  All  variables  in  (5.20)  save  ^ ' 

(2)  (2) 

and  satisfy  this  latter  requirement  trivially,  since  they  are  not 

time-differentiated  quantities,  and  d^^^  is  known^^  to  be  objective  and 

^kf^  readily  be  shown  to  be  objective  vectors  and  tensors, 

respectively.  To  see  this  consider 

where  Qj^£(t)  represents  an  arbitrary  time-dependent  proper  orthogonal  trans- 
formation and  arbitrary  time-dependent  translation.  In  (5,21)  tlie 

starred  quantities  represent  either  the  motion  as  seen  from  an  orthogonal  co- 
ordinate system  in  arbitrary  rigid  motion  with  respect  to  ours  or  the  motion 
plus  a superposed  rigid  motion  as  seen  from  our  coordinate  system.  From  (5,21), 


we  obtain 


(5.22) 


the  material  time  derivative  of  which  yields 


. Q^X^’/dt  + (de^,/dt)„<“ 


(5.23) 


Taking  the  spatial  gradient  of  (5.23)  and  employing  (5.21)^,  we  find 

5 (d»^"  */dt)dy;  . (dwf  >/dt)/dy_^  + 

35 

Now^  employing  the  well-known  relation 


■ Cki^'jei 


24 


in  (5.23)  and  (5.24),  respectively,  and  using  (5.22),  we  obtain 

,(1)*,..  ..(1) 


,(D* 
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dw.'-'  /dt  - wr'  0).^  - ^ 


^ (dw^^^Vdt)  - 4 - 

dy  by 

m m 


v( 


af-  ’ 

■‘z 


(5.26) 


which  shows  that  ^ij^  ^ objective  vector  and  tensor,  respectively. 

Obviously,  in  the  same  way  we  may  readily  show  that  and  constitute  an 

objective  vector  and  tensor,  respectively.  Now,  the  quantities  on  the  left-hand 

sides  of  (5.20)  cannot  be  arbitrary  functions  of  the  variables  shown  because 

arbitrary  functions  of  the  variables,  shown  will  not  satisfy  the  principle  of 
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material  objectivity  ’ , which  requires  the  constitutive  equations  to  transform 

appropriately  under  proper  orthogonal  treinsformations . However,  if 

D„(l)  D^(2)  D„(2)  , • ..u  ir 

<3 . , D.  . , tf.  emd  q.  are  expressed  in  the  forms 

J ij  J 1 


DtS 


- m D„(l)  a(1)  D^:(1)  *(1) 

ij  ^i,K?^j,L^KL  ' °ij  “ ^i,K^j,L^KL  •'  “^j,K^K  ^ 


D„(2)  .(2)  a,(2)  ,(2) 

°ij  " ^i,I^j,L^KL  ' j ~^j,K^K  " ‘^i“^i,K^K’ 


(5.27) 


where  T , and  L,,  are  functions  of  the  variables  shown 

on  the  respective  right-hand  sides  of  (5.20),  it  may  readily  be  shown  using 
established  methods^^  that  the  principle  of  material  objectivity  is  satisfied 
if  2nd  L are  vector  invariants  in  a rigid  motion  and  T and 

A^^  are  tensor  invariants  in  a rigid  motion.  Then  the  theory^^  of  invariant 

KLi 
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functions  of  vectors  and  second  rank  tensors  shows  that  the  required  invariance 

,(1) 

KL’  KL  ^ 'K  ' 'KL  ''  *K  K' 


is  assured  if  T._  , A^^,  respectively,  are  of  the  form 
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where 


(5.28) 


.(1) 


a)  .(1) 
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Sn  • ^i,M^j,N^ij  “ ^m  “ ’ S ’ 

2^^^  ■ V V B^^^  •=  y G = 0 (5  29) 

^i,M^j,N^ij  ’ “m  ^ijM^i  ’ S °,M  ’ ‘ ' 


^MN’  ^MN^  defined  in  (5.12)  and  (5.10).  Now^  it 


MN  ' M 


must  be  remembered  that  although  the  dependence  of  T_  , 

KXi  jQj  ' KXj  K 

and  L„  on  E.„.,  . P.^\  and  p/^,^  is  arbitrary,  there  are  conditions  on 

K MN’  M ’ MN  ’ M MN  ^ ’ 

tlieir  dependence  on  Zj^\  and  on  account  of  the 

Clausius-Duhem  inequality  (4.17).  Thus  the  dissipative  constitutive  equations 

in  the  general  case  are  given  by  (5.27)^  with  (5.26). 

Equations  (4,11),  (4.8)2_5,  (5 .13)  - (5 .15) , (5.27)  and  (5.28)  determine 

the  constitutive  equations  for  our  combined  continuum.  Thus,  all  that  remains 

in  the  determination  of  explicit  constitutive  equations  is  the  selection  of 

specific  forms  for  i|f,  T , A^^ , A^^ , and  L . Once  the  constitutive 

KL  K K K 

equations  have  been  determined,  we  have  a determinate  theory,  which  by  appropriate 
substitution  can  readily  be  reduced  to  10  equations  in  the  10  dependent  variables 
y^,  equations  are  the  three  each  of  (3,16),  (3.20), 

(3.25)  and  (4.16).  In  order  to  have  a complete  field  theory,  the  boundary  (or 
jump)  conditions  at  moving  surfaces  of  discontinuity  have  to  be  adjoined  to  tlie 
aforementioned  system  of  equations.  These  boundary  conditions  are  determined 


in  the  next  section. 


26. 


B^jfore  ending  this  section  it  is  perhaps  worth  noting  the  physical  fact 

that  the  objective  tensors  ^ij^^  nothing  more  than  the 

portions  of  dwf^Vdt,  (dwf^V^t)  dwf^Vdt  and  (dwf^Vdt)  , respectively, 

D j j 

beyond  that  of  each  due  to  the  local  rigid  body  rate  of  rotation  . As  a 

consequence,  if  the  aforementionec.  vectors  and  and  tensors  and 

■'  3 3 3,^ 

i ' 

w.  ! are  rigidly  fixed  in  the  continuum,  i.e,,  with  respect  to  then  the 
3 i ^ 

attendant  attenuation  will  vanish. 


6 . Boundary  Conditions 

In  this  section  we  determine  the  boundary  conditions  which  must  be  adjoined 

to  the  system  of  differential  equations,  as  noted  in  Section  5,  in  order  to 

formulate  boundary  value  problems.  As  usual,  these  boundary  (or  jump)  conditions 

are  determined  by  applying  the  integral  forms  of  ths  pertinent  field  equations 

to  appropriate  limiting  regions  surrounding  the  moving  (not  necessarily  material) 

39 

surface  of  dis'co.rtinuity  with  normal  velocity  u^,  and  assuming  that  certain 
variablei  rem»in  tyiunded.  The  pertinent  integral  forms  are  (3,2),  the  integral 
fomt  cf  (3.20),  (3.25)  eind  (4.17),  which  take  the  respective  forms 

(6.1) 
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pT]  dV  + J = J or  dV  ^ 0 , 


(6.2) 


(6.3) 


(6.4) 


where  pf  is  defined  in  (4,17). 
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For  all  the  integral  forms  considered,  a volumetric  region  is  taken  in  the 
39 

usual  waj,’  , and  it  is  assumed  that  all  pertinent  variables  remain  bounded . 

The  jump  conditions  obtained  from  the  respective  integral  forms  consisting  of 
(3.2)  and  (6.1)-  (6.4)  are 


u [p]  - n.  [v  p)  = 0 


n.  [T  . .]  + u [pv,]  - n.  [v.pv  ,]  =0 
i~  13~  n~  i~  1 ;j~ 


dll 


(1; 
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riL  dt  J 

Li  P 

dt 

(1) 
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(6.5) 

(6.6) 

(6.7) 

(6.8) 


n.[q./6]  - u [pTll  + n.fv.ptl]  SO  , (6.9) 

wh'^re  v'tt  have  introduced  the  conventional  notation  [C.)  for  cT  - C.  and  n.  de- 

X "^  XI  X 

notes  the  comp:nents  of  the  unit  normal  directed  from  the  - to  + side  of  the 
surface  of  discontinuity.  If  the  surface  of  discontinuity  is  material 

u = n.v^  - n.v.  , (6,10) 

n 11  11^ 

then  (6.5)  evaporates  and  (6,6)-  (6.9),  respectively,  reduce  to 


n.  =0  , 

n.  =0  , 

i~  ij  ~ ' 


u.  (o./ei  = 0 , 
1 "1 


(6.11) 

(G.)2) 

(6.13) 

(6.14) 


Moreover,  if  0 is  continuous,  i.e., 

[0]  = 0 , (6.15) 

across  the  .surface  of  discontinuity,  T is  bounded  and,  from  (4.17),  in  place  of 
(6.14),  we  have 
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n.  [q.]  "«  0 . 


(6.16) 


Wiis  latter  situation,  consisting  of  the  jump  conditions  (6.11)-  (6.13),  and 
(6.15)  - (6.16),  is  the  most  common,  and  if  the  body  does  not  abut  another  solid 
body  but  abuts,  say,  air  instead,  the  boundary  conditions  are  fully  defined  by 
tiie  noted  equations.  However,  if  a body  does  abut  another  solid  body  and  the 
full  field  equations  have  to  be  satisfied  in  each  region,  additional  conditions 
on  [y] , and  have  to  be  satisfied  at  the  surface  of  discontinuity. 


The  conditions  are  usually 


[yi  = 0 , = 0 , = 0 , 


(6.17) 


the  latter  two  of  which  may, by  virtue  of  (3.24)  and  (3.26)^,  respectively,  be 
written  in  the  form 

[(It  r<^^w<^U  r<2)w<^>]  =0,  ((It  r<2>)w<=>  . r<^^w(^)].0.  (6.18) 

«*w  ' r\^  1^/ 


Frequently,  the  thermal  conditions  are  such  that  v.u  may  eliminate  either  (6.15) 
or  (6.16).  Clearly,  all  boundary  expressions,  which  are  not  prescribed,  may  be 
expressed  in  terms  of  the  same  10  field  variables  as  the  10  equations  mentioned 
at  the  end  of  Section  6 by  making  the  appropriate  straight-forward  substitutions , 
We  can  determine  an  energetic  jump  condition  from  (4.4),  which , although 
not  needed  in  the  solution  of  many  types  of  boundary-value  problem,  can  be 
useful  for  obtaining  certain  types  of  information.  This  jump  condition  is 
obtained  by  applying  (4.4)  to  the  aforementioned  volumetric  region  surrounding 
the  (not  necessarily  material)  surface  of  discontinuity  and  assuming  that  all 
pertinent  variables  remain  bounded,  with  the  result 


dw. 


(1) 


dw, 


ji_  jk  k jk  dt  jk  dt 


(2) 

\ q.  I + u [T+  pe]  - n.  (v.  (T+  pe)]  = 0 , (6  .19) 

t ij  j ~ 
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where  T is  given  in  (4,3).  if  the  surface  is  material  we  have  (6,10),  and 
(6.19)  reduces  to 

dw^2) 

nfr  V (6.20) 

jL  jk  k 3k  dt  3k  dt  3_J 


7 . Generalization  to  N-Constituents 

In  this  section  we  generalize  the  equations  which  have  Ixicn  derived  for 
the  three-constituent  composite  material  to  a composite  with  N-constituents . 
Since  the  form  in  which  the  equations  for  the  three  constituent  composite 
material  have  been  written  makes  the  form  of  the  equations  for  th.-  N-constituent 
composite  rather  obvious,  we  briefly  present  the  basic  and  essential  resulting 
equations  here  for  completeness  without  presenting  any  of  the  intermediate 
equations.  In  faet,  where  possible  we  sir^ly  refer  to  the  generalization  of 
existing  equations  without  writing  new  ones . 

The  equations  of  Section  2 remain  unchanged  except  for  (2.3),  (2,9)  and 
(2.11)-  (2.13),  v/hich,  witli  the  exception  of  the  intermediate  equation  (2,11), 


take  the  resulting  forms 
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P„-  f,  pf’,  P - p'“>,  y pW„<">.o,  p<">  aw'»>/at-o.  (7,1) 

m=l  m=l  iii=l  TT.=1 

The  equations  of  the  conservation  of  linear  momentum  for  each  constituent,  i.e., 

(3.10)-  (3,12),  remain  the  same  except  that  the  sums  of  the  internal  interactions 

in  each  equation  increase  in  number  to  (N-1) , where  the  meaning  of  is 

obvious  from  the  discussion  in  the  last  paragraph  of  Section  2 and,  of  course, 

the  number  of  such  equations  increases  to  N.  The  equation  of  the  conservation 

of  linear  momentum  for  tiie  combined  continuum,  (3.16),  remains  the  same  and  the 

difference  equations  of  linear  momentum,  (3.20)  and  (3.25),  take  the  form 
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^ dt  J 


30. 


where 
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The  equations  of  the  conservation  of  angular  momentum  for  each  constituent , 

(3.13)  - (3,15),  remain  the  same  except  tliat  the  sums  of  the  internal  interaction 
couples  increase  in  number  to  (N-1),  where  the  meaning  of  the  is 

clear  from  the  discussion  in  the  last  paragraph  of  Section  2 and,  of  course, 
the  number  of  such  equations  increases  to  N.  From  Eg. (3.27)  it  is  clear  that 
the  equation  of  the  conservation  of  angular  momentum  for  the  combined  continuum 
now  takes  the  form 
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The  equation  of  the  conservation  of  energy  (4.1)  and  the  definition  of  the 

kinetic  energy  density  (4.2)  remain  the  same  except  that  tlie  sums  increase  to  N, 

Equations  (4.3)  and  (4.4)  retain  the  same  form,  but  the  sums  increase  to  (N-1) 

and  all  possible  quadratic  mixed  products  occur  for  (N-1)  terms.  The  equivalent 

of  Eqs,(4.5)  and  (4,6)  take  the  form 

N 


t 


= V 

L ~ 

n=l 


n 


(7,5) 


d(n)  ^ ^(n)  ^ ^(n)^  (N)  ^ » n • D 

^ rw  rv' 


n = 1,2, 


(N-1).  7.6) 


Equations  (4.8),  (4.9),  (4,11),  (4.12)  and  (4,14)  remain  unchanged.  The  number 
of  and  in  Eqs.(4,8)  increase  to  (N-1),  The  sums  in  the  remaining 

equations  in  Section  4,  i.e.,  (4.7),  (4,10),  (4,13)  and  (4,15)-  (4.17)  increase 
to  (N-1)  . Trie  nuutiber  of  occurring  in  the  generalization  of  the 

equations  in  Section  5 is  (N-1) . Consecjuently,  the  generalization  of  all 
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equations  in  Section  5 is  obvious,  i.e.,  wherever  the  w^‘  occur,  there  are  (N-1) 
of  them  in  place  of  2,  including  all  conjugate  quantities.  Accordingly,  the 
definitions  occurring  in  (5,10),  (5.13)  and  (5.29)  must  be  increased  thus 





» (dw^'^Vdt)  ’ «dwf"Vdt-w<’'^u,  , n«l,2,  ....  (N-1) 

Kg  ] ,K  ij  j J J- 


V , n-1,2,  (N-1) 


(7.8) 

(7.9) 


and  the  const! Vutive  relations  occurring  in  (5.14),  (5.15)  and  (5.27)  must  be 
increased  tlms 

Mf  - ■-  ■ "•*“> 

■ S’!  K''i  A"'  , ■'-1.2.  ••••  <“-i)  . (’-11) 

and  the  associated  dependence  in  (5.11)  and  (5.28)  must  be  increased  to  suit, 
i.e.,  must  contain  the  variables  with  superscripts  from  1 to  (N-1). 

Tlie  jump  conditions  (6.5),  (6.6),  (6.9),  (6.11),  (6.14)-  (6.16)  and  (6.17), 
remain  unchanged  and  tiie  sums  in  (6.19)  and  (6,20)  increase  to  (N-1).  The 
remaining  jump  conditions,  i.e.,  (6.7),  (6.8),  (6.12),  (6.13),  (6.17)2_2 
(6.18)  take  the  respective  forms 


n.(D.^"^]  + u [p^"^dTlf''Vdt]  - n Iv  p^"^dT]f"Vdt]  = 0, 
1 j ^ n~  j ~ 1 3 

n,  = 0 , 

ig  ~ 

= 0 , 

^ <V 

N-1 

[wf">  t =0. 

~ g /_  g -V 

m=l 


(7.12) 


(7.13) 


(7.24) 


(7.15) 
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8.  Piola-Kirchhoff  Form  of  the  Equations 

Up  to  this  point  all  the  equations  have  been  written  in  terms  of  present 
(or  spatial)  coordinates.  Since  the  reference  (or  material)  coordinates  of 
material  points  are  known  while  the  present  (or  spatial)  coordinates  are  not, 
it  is  advantageous  to  have  the  equations  written  in  teniis  of  the  reference 
coordinates.  To  tliis  end,  analogous  to  the  Piola-Kirchhoff  stress  tensor  j ? 
which  is  defined  by 

n.T. . dS  « N K. . dS  , (8.1) 

1 13  X Lg  o ^ 

we  define  the  reference  relative  stress  tensors  by 


as  . as  , 

1 X Lg  o’ 


(8.2) 


Where  dS^  and  denote  the  magnitude  of  and  unit  normal  to  an  element  of  area 

in  the  reference  configuration,  which  has  magnitude  dS  and  unit  normal  n^  in  the 
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present  configuration.  By  virtue  of  tlie  well-known  relation 


n.  dS  = JX  .N  dS  , 

1 X,i  X o’ 

from  (8.1)  and  (8,2),  in  the  usual  way,  we  find 

K.  , = JX  ,T.  . , •>  JX 

Xj  L,i  Lg  L,i  13 


(8.3) 


(8.4) 


Now,  using  (8.4)^  and  (8,4)2,  respectively,  with  (3.16)  and  (7.2)  and  employing 
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(2.8).  (2.10)  and  tlie  well-known  identity 


Ve  obtain,  respectively, 


K . - + p f . = p dv./dt  , 
X3,l,  3 ’^o  3 ’ 


(8,5) 


(8,6) 


(n)~(n)  (n)  _ (n)  2 (n)  2 

Lj,L  Po  ‘j  + - Po  , n-1,2, 


. . . (N-1) 


(8,7) 
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the  first  of  vhich  is  the  Piola-Kirchhof f form  of  the  stress  equations  of  motion 
and  the  second  of  which  is  the  reference,  form  of  the  (N-1)  relative  stress 
equations  of  motion. 

Analogous  to  the  foregoing,  we  now  define  the  reference  heat  flux 


vector  Q by 
L 


(8.8) 


"A  • 

which  with  (8.3)  yields 

Substituting  from  (8.9),  (5,27),  (5.29)  and  (7.8)-  (7.10)  into  the  generalization 


of  (4.16)  for  N-constituents  and  employing  (8.5)  and  (2.10),  we  obtain 

(N-1) 

(m)„(m)  ,^(m)„(m) 

‘KL““KL^“‘'  L ' KL 
m=l 


♦ I - fiL.L-Po®  '“-I"’ 


Which  is  the  reference  form  of  the  dissipation  equation. 

In  view  of  (0.4),  (8.9),  (5 .13)  - (5 .15) , (2.10),  (5.27),  (7.10)  and  (7.11), 
the  pertinent  constitutive  equations  for  this  section  may  be  written  in  the  form 


(N-1) 

9'^  . 

V 

r 9't' 

pQ  L 

L.  (m) 

IjM 

m=l 

dN ' 

J-l 

0\j[ 

H r.\n) 

-^1.3 

SP. 


LM 


di. 


^o"j,K  (n)  ^ 
KL 


(8.11) 


= Jy  , Q = JL  , 

Lg  -^j,M  L.M  ' Y.’ 


(8.12) 


where 


K.  . = V . + V . , . (8.13) 

L3  L3  Lj  ' g g 3 ^ Lg  Lg  Lg  - ' 


I'U 


I 
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and  * T and  L arc  the  generalized  versions  of  (5.11)  and  (5.28), 

LM'’  K ■’  I.M  K 

which  are  discussed  in  Section  7.  In  view  of  (8.1),  (8,2)  and  (8.6)  - (8.8),  in 
reference  coordinates  tJie  boundary  conditions  (6.6),  (7.12),  (6.9).  (6.11), 
(7.13),  (6.14)  and  (6.16)  take  the  respective  forms 


[K.  .]  + U„P  [v  ] = 0 , 
Lj~  N o~  j~ 

(8.14) 

+ U„p  [dflf^Vdt]  =0,  n = 1,2, 
NO  9 

, . (N-l)  , 

(8.15) 

- Vol'Hl  • ° ' 

(8.16) 

(N-1)  , 

(8.17) 

(8.18) 

Where  U is  the  intrinsic  velocity  of  the  singular  surface,  i.e.,  the  velocity 
N 

of  the  singular  surface  in  the  reference  coordinate  system.  The  boundary  condi- 
tions (6.15),  (6.17)^,  (7.14)  and  (7.15)  remain  unchanged,  while  (6.5)  degenerates 
to  nothing  in  the  reference  coordinate  description. 


9,  Linear  Equations  for  the  Two-Constituent  Composite 

In  this  section  we  obtain  the  linear  eguat’Oi;s  for  the  two-constituent 
composite  material  in  the  absence  of  dissipation  fron  the  general  nonlinear 
equations  for  the  N-constituent  composite  material  j resented  in  Sections  7 and  8. 
To  this  end  we  first  note  that  for  the  two-constituent  composite  N = 2.  Then, 
in  the  usual  way,  we  define  the  mechanical  displacement  vector  Uj^  by 


(9.1) 


where  6,  is  a translation  operator,  which  serves  to  translate  a vector  from 
iM 

tlae  present  to  the  reference  position  and  vice-versa  and  is  required  for  nota 
tional  consistency  and  clarity  because  of  t)ie  use  of  capital  and  lower  case 
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indices,  respectively,  to  refer  to  the  referenct  and  present  coordinates  of  the 
center  of  mass  of  material  points.  From  (9.1),  wc  have 


^ 6.^  + 
iii 


6 u 
iM  M,L’ 


(9.2) 


and  substituting  from  (9.2)  into  (5.12)  and  neglecting  products  of  u 

M,  tj 

obtain 


E <=a  e 
LM  LM 


= 1 


we 

(9.3) 


which  is  the  usual  infinitesimal  strain  tensor.  .Similarly,  substituting  from 
(9.2)  into  (7.7)  for  N = 2,  we  obtain 


^(1)  o (1)  (1)  „(1)  . '1)  (1) 

P ssfi  w N W — w'^ 

LM  kL  k,M  L,  L k L ^ 


(9.4) 


where  we  have  taken  the  liberty  of  utilizing  capital  indices  to  denote  tlie 

Cartesian  components  of  tlie  relative  displacement  vector  w'  ' in  the  linear 

\IL.(3.'  R il) 

' and  are  assumed  to 

vanish  when  u and  w^^^  vanish,  in  this  linear  theory  ijr  must  be  a homogeneous 
M L 

quadratic  function  of  the  form 


,1  1 (1)  (1)  1 , (1)  (1) 
pt=—  c ee  + — aw  w +—  b v w + 
2 KLMN  KL  MN  2 KL  K L 2 KLMN  K,  L M,  N 


n \ 


ni 


a e w'*'  + B e w'  ' + y w'  'w"'  . 
MKt  KL  M '^KLMN  KL  M,  N 'mKL  K,  L M ’ 


(9.b) 


where  the  the  usual  elastic  constants  of  ordinary  linear  elasticity, 

the  a may  be  called  the  difference  displacement  elastic  constants,  the  b, 

lUj  ' KLMN^ 

the  relative  elastic  constants  and  ^KLMN  ^MKL’  respective  coupling 

constants.  In  the  arbitrarily  anisotropic  case,  tliere  are  the  usual  21  inde- 
pendent = 6 independent  , 45  independent  b_  18  independent  . 

54  independent  3^,..  and  27  independent  Y„,_  > for  a total  of  171  independent 


material  constants. 
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Now^  substititing  from  (7.3)^,  (8. IX),  (8,13)  and  (9.2)  ~ (9.5)  into  (8.6) 


and  (8.7)  and  neglecting  ail  nonlinear  terms,  we  obtaiis 


4.  -Sl 


1 + 

where  we  have  employed  the  relations 


(1)  ^o*^K  .M  Pq  M ’ 


V,  = Sy./5t  * 6 . 9u^St  = 6.  u , f:  = 6.„f„ 

2 3 ]M  K gM  3 ]M  M 


ad) 

J 3^!^^  = { 
3 

3(Pq'1') 

,(1) 

-.d) 

o 

CL 

1 

1 

II 

“ ■ 

From  (9.5)  and  (9.9)  ve  have 

“ ‘^LMKS®KN  '^KLM^K  ^LMKn'^K,  N > 

*^LM  ^KIIML^KN  '^RML^'^K  ^MLKn'^K,  N ’ 

*^M  ‘^MKL^KL,  ” ' MKl'^K,  L ’ (9-10) 

The  substitution  of  (9.10),  witli  (9.3),  into  (9.6)  and  (9.7),  respectively, 
yields 

'^LMKN^K,  NL  ’’’  '^KLm”k,L  ^ ^LMKn''V,  KL  ’ (9.11) 

^ u + y + b W^^^  + -— £ — p 

•^KNLM  K,NL  ’kLM  K,  L LMKN  K,  NL  , ^ (1)  '^O  M 

7 ; ^ 1 + r 

-a  u - a w^^^  - V w^^^  = r^^^p  w^^^  (9  12) 

“mKL  K,L  ^ML  L ^MKL  K,L  ^O  M ’ ( ■ > 

which  constitute  6 linear  differential  equations  in  the  6 dependent  variables 


u and  w'  . To  this  system  of  equation.?  we  must  adjoin  the  linear  boundary 
K K 
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conditions  across  material  surfaces  of  discontinuity^  which  are  obtained  by  sub- 
stituting from  (9.8)^  ^ into  (8.17)  for  N = 2 with  the  result. 

If  the  surface  abuts  space,  one  side  of  the  jump  brackets  in  (9.13)  determines 
applied  traction  terms  in  the  usual  way.  On  the  other  hand  if  the  body  abuts 
another  solid  body,  we  must  obtain  the  associated  jump  conditions  by  substituting 
from  (9.1)  and  (9,4)  into  (6.J7)^  and  (7.15),  respectively,  with  the  result 

[ (1  + = 0 . (9.14) 

The  linear  equations  for  the  two  constituent  composite  with  discontinuous 
reinforcement  (chopped  fiber)  Ccin  be  obtained  from  the  foregoing  equations  in 
this  section  simply  by  setting  b_„„,  and  Y_„  equal  to  zero  wherever  they 

jtJjMN'  KXirlN  KIiM 

occur.  Under  these  circumstances  vanishes  and  the  boundary  conditions 

LM 

{9.13)2  and  (9.14)  do  not  exist. 

10.  Material  Symmetry  Considerations 

In  this  section  we  obtain  the  linear  equations  for  the  isotropic  and  trans- 
versely isotropic  two-constituent  composite  material.  Although  we  can  obtain 
these  equations  directly  from  the  arbitrarily  anisotropic  equations  presented 
in  Section  9 by  writing  the  tensors  for  the  particular  symmetry  involved,  it  is 
advantageous  to  return  to  the  stored  energy  function  p^iti  and  write  all  the 
quadratic  scalar  invariants  first  for  the  isotropic  material  and  then  for  the 

transversely  isotropic  material,  especially  when  a great  deal  of  symmetry  exists, 

43  44 

as  in  tliese  two  cases.  The  tables  of  integrity  bases  provided  by  Spencer  ’ 
for  the  two  transformation  groups  involved  prove  to  be  extremely  valuable  in 
obtaining  the  independent  quadratic  invariants,  we  are  concerned  with  the 
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quadratic  invarieint.s  of  a symmstric  tensor  an  asymmetric  tensor  and 

X»n  n 

a vector  Since  Spencer  systematically  considers  syircnetric  tensors, 

Li 

antisymmetric  tensors  and  the  skew- symmetric  tensors  of  vectors,  we  must  de- 

coiMxjse  the  asymmetric  tensor  into  its  symmetric  and  antisymmetric  parts 

I4,  M 

p^  and  , respectively,  and  write  the  skew-symmetric  tensor  of  the 

LM  LM  liM 

vector  thus 

K 


.(1)  . ..a) 


.a)  ..(!)- 


n »-  fw'  ' + w'  M D =—  <w'  ' -w' 

^LM  2 'l,M^  M,L^^  ^LM  2 ^ L,  M M,l'  ’ 


W = e w<^^ 
LM  LHK  K 


(10.1) 

(10.2) 


Then  can  be  written  as  a quadratic  polynomial  invariant  in  the  sum  of  the 

. . S 

invariant  products  of  two  symmetric  tensors  e and  ^ an  antisymmetric 

XiM  Xifr 

tensor  pf.  and  a skew-symmetric  tensor  W ... 

LM  LM 

For  isotropic  materials  possessing  a center  of  symmetry,  is  a scalar 
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Invariant  under  tlie  full  orthogonal  group.  Spencer  lists  the  basic  invariants 
of  a number  of  second  order  symmetric  and  antisymmetric  tensors  for  the  proper 
orthogonal  group.  Prom  this  list  all  quadratic  invariants  under  the  full 
orthogonal  group  may  readily  be  obtained.  Thus  we  find  that  for  an  isotropic 
material  the  homogeneous  quadratic  function  p^i|(  may  be  written  in  the  form 


A 

T e 


2 'KK^LL  ^I^KkPiJ,  ^2®KL^LK  2 ^1""KK"LL 


P P 

ir  iz-trC"  T 


1 

3^KL^LK  2 


A A 
+ b-,P„P. 


(1)  (1) 


(10.3) 


Substituting  from  (10.3)  into  (9.9)  and  employing  (9.3)  and  (10. L),  we  obtain 


39. 


^LM  “ ^l^K,  K®LM  '*'  2 ^2^'^,M  '''  '^M.L^  "*  K®LM  "''  *^2  m'*' ''m,  L^ 


+ b,  (w, 


(1)  ..(1) 


3'  L 


' - w'  M 
, M M,l'  ’ 


- - a w<^> 
M 1 M 


(10.4) 


which  are  the  linear  constitutive  equations  for  the  isotropic  two-constituent 
composite  material.  Substituting  from  (10.5)  into  (9.6)  and  (9.7)  and  ignoring 
the  exLernal  body  forces,  we  obtain 


KK  * (^1  ^ 2 ^2Vk,  KM  2 KK  *^o^M  ’ (10.5) 


V^l  '*'1^2)  ■*■  2 ^2%  KK"*"  ^1'^^2'^^3^'^K, 


(1) 
KM 


'2  "3'  M,  KK  1 M 


(10.6) 


which  are  the  equations  of  motion  for  the  isotropic  two  constituent  composite 
material. 

A material  with  a single  preferred  direction  which  is  the  same  at  every 

point  is  said  to  be  transversely  isotropic.  For  such  a material  p^|  is  a 

scalcur  invariant  under  rotations  about  the  preferred  direction,  which  we  take 

along  the  x^-axis.  The  transformations  under  which  we  liave  invariance  are  the 

rotations  about  the  x^-axis,  and  reflections  in  the  planes  containing  the 
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axis.  Spencer  gives  a list  of  invariants  under  this  transformation.  In 

addition  to  the  foregoing  we  require  invariance  under  reflections  in  the  plane 

44 

perpendicular  to  the  x^-tucis.  From  Spencer's  list  all  quadratic  invariants 
satisfying  the  latter  additional  symmetry  requirement  may  readily  be  obtained. 
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Ihus  we  find  that  for  our  transversely  isotropic  material  the  homogeneous 
quadra?;ic  function  can  be  written  in  the  form 

'’o*  ” 5 'I'ao'es  * °2'm‘33  * ■"  =4'3t.'<33  * Y 's'33'33 

" ■"  ^3^33  * ^4*30^^  * ^'33^^ 


« S j'SS  .Ir^SS 

h^3-^^32  ^ ^ ^ '^2Wpe  ^ V<«^33  ’^4p3oPa3 


s s 


2 “5^33^33  ^7^3c^3o  ■*■  2 l a a 2 2 3 3 , 

(10.7) 


where  the  greek  indices  take  the  values  1 and  2 and  skip  3. 

Substituting  from  (10.7)  into  (9.9)  and  employing  (9.3)  and  (10.1),  we 
obtain 

y^ap  * °2'^3,3*ap  ^l^Y,  W 

^3  3,3  a'3  2 P2'  a,P  p.a  ’ 

K - = K_  = C.  (u.  +u  .,)  + -:r  P.,  (w'  ^ + w'  ,)  , 
a3  3a  4 3,  a a,  3'  2 '^4  3 a a,  3 ’ 


K = C U + ^r'^-3  -3  ■*■  ■*■ 

.33  2 a,  a 5 3,3  5 a, a 6 3,3’ 

' Vv,-(  ‘aP  *1  ^2%,p  * '‘p,.’  * \"3,3*»p  + 


.(1) 


.(1) 


.(1), 


.(1>  ..(1) 


^ v;?ap  ^ V'3,?ap  " - %P>  ’ 

tu  .)  tb,(wf>  .w<’-^  .ib  rwf>-w<^>), 
a3  t 4 3,a  a,  J 4 3,  a a,  3 2 7 3,  a c^,  3 


v(l)  1 


.(1)  . ..(^).  . 1 c:  ,..(1)  ..(1) 


= - p4  (u,  „.  + u^,  ,)  + b,  (W3^^;  (W-'  - W3- ) , 


'3a  2 ^^4'  3,a  a,  3'  "4'  3. 


(i) 


t b_w. 


(1) 


■^33  " ^3% a ^ V’3,3  ^ “5"3,3  ’ 


..(1)  ' (1) 

^ p ■ “ ^I'^y  ’ 


.(1) 


A (1) 

®2’"3  ’ 


(10.8) 
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which  are  the  linear  constitutive  equations  for  the  transversely  isotropic  two- 
constituent  con^site  material  witii  the  preferred  direction.  Substituting 
from  (10.8)  into  (9.6)  and  (9.7)  and  ignoring  the  external  body  forces,  we  obtain 


^l'^a,Qfp  3p  Pa  ^ ”p,oa^  '^^'^ajaP  ^ 2 pa  ^ '^PjOa^ 

®3''3,  3p  ‘"4^'^3,P3  '^p,33^  '*'2  ^4^'*3,P3  ''p,  33^  “ "^O^P  ’ 


« - , , 1 a , (1)  (1)  V u (1) 

c„u  ^ + c,  (u,  + u n ) + T P.,  + w'  ) + c-  U-  ,-  + Pc'^  , 

2 a,a3  4 3 oa  a,  3a  2 4 3,aa  a,  3a  5 3.33  5 a o3 


•*  V3^33  - ’ 

Va,ap  2 "^2%aP  * "p,oa’  2 ^4<"3,3p  '"p,33) 

^5^3P  " ^<ap  " "P,aa^  " Va,aP  " Vi,3P  " ^e^L'^Pa 


(10.9) 


P3  ■*■  '*p,33^  "^2  ^7^''P,  33  ■■  '’3,  3p^  “ ^iS  P 


P ' 


(1)  (1)  ..(1) 


o P 


§_u  , + 4 P.  (u_  + u _)  + P,u-  + b.,w'^\  + b.  (w'"*^^  +w'‘*’^,) 

*^3  a, a3  2 ^4'  3,0fa  a,a3  ^6  3,33  3 a, a3  4 3,00  o, o3 


(1) 


.(1) 


. 1 CL  , (1)  (Da  d (D  ^ (D  (D  -(D 

+ -g  b„  w'  - w _ ) + b,-W_  - a_w  * r p w 

2 7 3,00  a,  3a  5 3,33  2 3 *^o  3 


(10.10) 


which  are  the  equations  of  motion  for  the  transversely  isotropic  two-constituent 
composite  material  with  x^  the  preferred  direction.  We  have  boulierod  to  write 
the  linear  constitutive  and  differential  equations  for  the  transversely  iso- 
tropic two-constituent  con^josite  material  in  complete  detail  because  we  deem 
this  to  be  a particularly  important  symmetry  for  fiber  reinforced  composites. 


42. 


11.  Wave  Propagation 

The  solution  for  plane  wave  propagation  in  the  arbitrarily  anisotropic 
infinite  medium  may  readily  be  obtained  by  substituting  plane  waves  in  the  system 
of  linear  equations  in  (9.11)  and  (9.12).  However,  the  resulting  algebra  is 
sufficiently  lengthy  and  cumbersome  that  it  becomes  relatively  involved  to 
extract  useful  physical  information  from  the  resulting  system.  Since  the  iso- 
tropic case  contains  meiny  of  the  interesting  features  concerning  the  propaga- 
tion of  waves  iii  the  two-constituent  composite  and  is  much  less  cumbersome 
than  the  general  anisotropic  case  because  of  the  considerably  smaller  number  of 
material  constants,  we  treat  plcine  wave  propagation  in  two-constituent  isotropic 
composites  in  this  section.  To  this  end  as  a solution  of  (10.5)  and  (10.6) 

consider  ^ 

i ( Sn.  A -iiit ; i V ■pu  X -Ci}L > 

u.  = A.e  , w.  - B.e  ^ , (11.1) 

3 2 ’33  ’ 

where  is  a unit  vector  denoting  the  wave  normal.  The  solution  (11.1) 
satisfies  (10.5)  and  (10.6),  provided 

(p  §2  _ ^ (X  +^,)  ^\n^n.  + I 

+ (01  + = 0 ,1 
^ ^ a.  j 

I * (Pi  + I P2>5\Vj*[<‘’2-*’3>5'  + »i  - -^iPo 


+ 0=1  +*>2  Wj  ‘O' 


(11.2) 


At  this  point  it  should  be  noted  that  in  order  to  secure  the  positive  definite- 
ness of  p^i|r  in  (10.3),  we  must  have  the  conditions 

p>0,  3X+2|J.>0,  b.,>0,  3b^  + > 0 , b^  < C , > 0 . (11.3) 

Equations  (11.2)  constitute  six  linear  homogeneous  algebraic  equations  in  the 
A.  and  B.,  which  may  also  be  regarded  as  two  vectorial  equations  of  the  form 
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■ 0 aniS  S',  “ 0.  since  the  medium  is  isotropic,  a major  reduction  in  the  algebra 
results  if  v;e  decompose  each  vector  equation  along  the  direction  of  the  unit 
normal  n and  in  the  plane  normal  to  n.  To  this  end  we  write  A and  B in  the  form 


A*»n  °An  - nXnXA,  B*n*Bn  - nXn 

•’V  <V/  ^ 


X B 


(11.4) 


and  first  we  write  the  equations  in  the  normal  direction  n,  to  obtain 

((X  + 2^)5^  - p^(U^]n  • A + (0  +fi  )5^n  • B = 0, 

2 2 2 
(8,  + g_)§  n • A + t(b,  +2b.)§  +a  - r'  ’pOi  ]n  • B = 0,  (11.5) 

which  constitute  a system  of  two  linear  homogeneous  algebraic  equations  in 
n • A and  n • B,  tlius  showing  that  purely  lonaitudinal  waves  exist  in  tlie  iso- 
tropic  two-constituent  composite.  For  a nontrivial  solution,  the  determinant 
of  the  coefficients  oi  n ' A and  n • B must  vanish,  which  yields 


.(1)  2 4 


u.'‘‘  - f[(bj^  + 2b2)p^  + (X  + 2^i)  + 

+ t(<'^  + 2|i)(bj^  + 2b2)  - + 2p,)aj^]5^  = 0 


(11.6) 


Equation  (11.6)  governs  the  propagation  of  longitudinal  waves  in  the  isotropic 

two-constituent  con-iposite.  Prom  (11.6)  it  is  clear  that  for  a given  wavenumber  5 
o 2 

there  are  two  ou  , and  solving  tor  U)  and  expamding  for  small  we  obtain 


ai^  - “(X  + 2p. ) ^^  + 0 ( ) , ui"  = —7^ 
1 Po  2 „(1) 


2 ^1  ^’^l  ■*'  ^^2^  ^2  . _ .,4, 


. . (1) 
r p r p 
o o 


5 +0(g  ) . 


(11.7) 


In  view  of  (11.3),  we  see  from  (11.7)  that  both  a-^  and  are  real  for  real  |. 

Moreover,  it  is  clear  that  on  an  u)  vs  g diagram  tliere  are  two  branches,  one 
emanating  from  (U  = 0,  g - 0,  with  positive  initial  slope  ^(X  + 2p.)/pJ,  and  the 
other  emanating  from  O)  = “ 0,  with  zero  initial  slope  and  positive 

curvature  (b^  + 2b^ ) /\^  • 

Now  we  write  tlie  equations  in  the  plane  normal  to  n to  obtain 
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^ ~ p w^)nXnXA  + -|  P^l^nXnXB  = 0,  i fi  §^nXnXA  + [ (b  - b ) 5 ^ 

Q rw  r\d  r\t  ^ ^ ^ ' J J J 4 


2 '^2- 


2 3' 


(1)  2 

+ a,  - X puj]nXnXB«0, 
X O ^ 


(11.8) 


which  constitute  a system  of  two  linear  homogeneous  algebraic  ec[uations  in  tlie 
vectors  n X n X A and  n X n X B,  thus  showing  tliat  purely  transverse  waves  exist 
in  the  isotropic  medium.  For  a nontrivial  solution,  the  determinant  of  the 
coefficients  of  n x n X A and  n X n X B must  vanish,  which  yields 


r^^^p^  UJ^  - [[r'^^^p  p,  + p (b  - b-)]  5""  + a p lu)'" 

O O 0^*3  XO 


(1) 


+ [tp  (b2~b3)  - i tpa^H^  = 0 . 


(11.9) 


Equation  (11.9)  governs  the  propagation  of  transverse  waves  in  the  isotropic 

2 

medium.  From  (11.9)  it  is  clear  that  for  a given  v.’avenumber  f there  are  two  ui'  , 
2 

and  solving  for  (U  and  expanding  for  small  5,  we  obtain 

^2  U 2 4 

\ o(r) 


.^2  _ ^^2  " ^3^  ,,2  ^ 

^2--(Tr^~(i) — ^ 

^ Po  ^ Po 


(11.10) 


Again,  and  tor  the  same  reasons,  it  is  clear  UiaL  both  and  are  real  for 

real  5;  an  cu  vs  ^ diagram  there  are  two  branches,  one  emanating  from 

d)  = 0,  § = 0,  with  positive  initial  slope  ,/li/p^,  and  the  other  emanating  from 
I nj 

UJ  = p^,  1 = 0,  with  zero  initial  slope  and  positive  curvature 


(b„ 


b3)//'^^Po^l  • 


12 . Dynamic  Potentials 

In  the  classical  theory  of  isotropic  linear  elasticity  it  is  possible  to 
reduce  the  displacement  equations  of  motion  to  wave  equations  in  the  Lame 
potentials  by  means  of  the  Helmholtz  resolution.  In  this  section  the  analogous 


reduction  is  obtained  for  Eqs. (10.5)  and  (10.6),  which  were  derived  for  an  iso- 
tropic two  constituent  composite  material,  and  the  con^leteness  of  tlie  repre- 
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sentation  is  e.stabllslied  following  Sternberg  with  minor  modifications.  To 
this  end  we  write  Eqs. (10.5)  and  (10.6)  in  the  invariant  vector  form 


(12.1) 


(p,u  + ~ p w)  + VV  ■ I (X  +p)u  + (p  + i P,)w]  = pii 


where  we  have  taken  the  liberty  of  omitting  the  superscript  (1)  from  w and  r. 
The  substitution  of  the  Helmholtz  resolutions 


u = V0,  + V X H-  , V 


w = vx,  + V X G , V • G 


0, 


(12,3) 


into  (12,1)  and  (12.2)  yields 
V[(X  + 2u)5\t 

V[(P^  + P^)V\+  (b^  + 2b,)»\-.^X^-rp^X^!  + 

J [i  ®2''^Sl  ► l»2  - »3 ' - *lSl  - *Po£l]  - “ ■ <12  ■ 5 > 


Taking  the  divergence  and  the  curl,  respectively,  of  both  (12.4)  and  (12.5), 
we  obtain 

(X  + 2^I ) + (P^  + - pj^]  = 0 , (PV^H^  + 1 p^v^G^  _ p = 0 , 

+ 2b2)v\^  - - rp^Xj^]  = 0 , 


''11  <h-l'3’'''Sl-hSl-'PoSl]-“- 


(12.6) 


where  we  have  employed  (12.3)  along  with  the  identity 


VXVXV  = VV*V  - VV, 


/Xy 


(12.7) 
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in  arriving  at  (12. 6).  Clearly,  from  (12.6),  we  may  write 

(l  + 2|r)V^0j  + i - p = b , 


where 


(P^  + ~ ^Po\  = ^ y 

i * (>'2  - »3>’"Sl  - nSl  - ‘Pjl  ■ J . 


V^a  = 0,  V^b  = 0,  Vb-0,  V^d » 0 , V^e  = 0 , V*e  = 0. 


(12.8) 


(12.9) 


Now,  let 


+ x = X3^  + d,  g = g^  + e, 


(12.10) 


where  A,  B,  D aind  E are  particular  functions  to  be  selected  in  order  that 


pA  = a,  pB  = b,  a,D  + rp  D = d,  a, E+rpE  = e 
o ’ cw  ~ ’ 1 o ’ 1~  o~  ~ 


(12.11) 


To  this  end  we  take  A,  B,  u and  E in  the  forms 


A a -i-  I a (X,  s ) dsdT  , B = -i~  [ [ b (X,  s ) dsdT  , 

P 0 J ' rsi  P V V ~ ^ 


D = d (X,  s)sin  o)  (t  - s)ds  , E = | e (X,  s)sin  UJ  (t-s)ds  , 

u)  rp  J o ' ~ oi  rp  J o ’ 

o o o 


I rp  J 
o o o 


(12,12) 


w " a /ro 
o 1 *^0  ’ 


(12.13) 


and  on  account  of  (12.9) 


V^A“0,  V^B  = 0,  V*B  = 0,  V^D  = 0,  V^E  = 0 , V*E  = 0.  (12.14) 

^ cw  ' <v,  ' r\j  * 2-u  rw 


Moreover,  from  (12.9)^_^  and  (12.11)^  ^ we  have 


V^D  = 0 , V^E  = 0 . 


(12.15) 


Substituting  from  (12.10)  into  (12.8)  and  employing  (12.11)  and  (12.14),  we  obtain 


1 


(X+2UL)V^e+  (3,  + P,)v\-p  0 = 0,  + i p V^G-p  H = 0. 

i.  4L  O Q~ 


(Pi  + ^2^  2b,)v\  - aiX  - rp^X  - 0 , 

i P-V^H  + (b^  - b,)V^G  - a ,G  - rp  G = 0 , 

2.  2.  2 2i  0~ 


From  (12.3)  and  (12,10)  we  may  write 


(12 . 16 ) 


where 


u=V9+VXH+u,  w = VX  + VXG+w, 


u = -VA-VXB,  w = -VU-VXE. 


(12.17) 


(12.18) 


Equations  (12.18)  and  (12,14),  along  with  the  identity  (12,7)  applied  to  B 


and  E,  imply  that 


V‘u=0,  Vxu=0,  7*w=0,  VXw=0. 


Therefore,  there  exist  functions  C(X,  t)  and  T](X,  t)  such  that 

u = VC  , V^C  = 0 , w = VT] , V^Tl  = 0 , 
and  (12.17)  and  (12.20)  permit  us  to  write 


u = ve+VC  + VXH,  w = VX+VTi  + VXG 


(12.19) 


(12.20) 


(12.21) 


]^]7oin  xii'cu  diici  anu.  cin^ioy^in^  (1-2,  j_w)  ctrivl. 


(12.20)^  we  obtain 


yC  = 0,  V(a /n  + rp  h)  = 0 

r*,i  ' ^ 


(12.22) 


Hence,  C ^*^‘3  I]  must  have  the  respective  forms 

C = Q'j^(t)  + t,P^(X)  + y^(X)  , 7]  = (t)  + P^  (X)  cos  uu^t  + (X)  sin  lu^t  . (12.23) 


Equations  (12.20)  and  (12.23)  indicate  that 

V^P^  = 0,  V^y^=0,  V^p2=0,  V^^(2  = 0 


(12. 24) 
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Finally,  let  us  define 


cp  = 0+  C-Oj_,  'j'=X+11-a2, 


(12.25] 


which,  with  (12.21),  (12. 3)^  (12.10)^  and  (12.14)^  ^ enables  us  to  write 


u = Vcp+VXH,  V*H  = 0, 


w = V\|(+'7XG,  V-G  = 0, 

f\f  f\f  /\j  /V<  ■’  r\j  /V  ' 


(12.26) 


and  from  (12,16)  and  (12.23)-  (12.25)  we  see  that  the  potentials  satisfy 
(\  + 2M.)V‘'cp  + (P  + 3 )V^t|;  - p 9 = 0 p,V^H  + ^ P-V^G  - p H = 0 , 


X ? 2 • • 

— p V H + (b„  - b ) V G - a G - rp  G = 0 , 
2 2 ~ 2 3 ~ 1~  o~  ^ 


(12.27) 


and  wc  have  shown  that  the  representation  is  complete,  it  should  be  noted 

given  by  (12. 26)^^  satisfy  (12.1)  and  (12.2)  provided 
(12.27)  hold  even  if  H and  G are  not  solenoidal,  as  may  be  confirmed  by  direct 
substitution. 


13 . static  Potentials 

A complete  solution  of  the  displacement  equations  of  equilibrium  for  an 
isotropic  tv/o-constituent  composite  material 


V^(p,u  + i P w)  + V7  • I (\  + ;j.  )u  + i 3,  + i 3.,'Vvl  + 0 f = 0 , 

\~  ^2'^/  ~~L  'N.  '^x  2 2/~J 


(13.1) 


- a,w  + [r/(l  + r)]  p f = 0 , 
i~  '^o~  ’ 


(13.2) 


is  obtained  in  terms  of  stress  functions,  which  reduce  to  tlie  Pajikovitcli 
functions  of  classical  elasticity . Ihe  procedure  follov;s  that  of  Mindlin"^^’ 


^ - ■(■si-’ 
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Tl\e  substitution  of  the  Helmholtz  resolutions 


u = Vcp  + VXH,  V*H  = 0, 


w=^Vil(  + VXV,  V*V  = 0, 


(13.3) 


into  (13.1)  and  (13,2)  yields 


V ta,Vcp  + V X H + x v)  + f = 0, 

+ + (b^  + 2b2)V^t  - a^f]  + 

+ f = f) , 


vxfi  + (b,,-b  )V‘V 

~ L2  2 ^ 2 3 ~ 


(13.4) 


(13.5) 


where  for  convenience  we  have  introduced  the  definitions 


f = p f , f [r/(l  + r)]p , 

~ O'-*  O'-  ' 

a = ^ ^ a = 


M. 


M.  ' 3 2p. 


(13.6) 

(13.7) 


Let  us  define  a vector  function  B by 


B = a,Vcp  + V X H + U'-Vt  + Q\,V  X V, 


then 


Taking  the  divergence  of  (13,8)^  we  obtain 


(13.8) 


(13.9) 


• B = V (Of^cp  + Qf-'f)  , 
~ X ^ 


(13.10) 


which^  with  the  definition 


X = 0^9  + »2'1'  , 


(13.11) 


enables  us  to  write 


V ■ B = V X . 


(13.12) 
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Since 


1 2 
V (r  * B)  ■-=  2V  • B + r * V B , 

r\j  «S#  ' 


(13.13) 


from  (13.12)  and  (13.9)  we  have 


V B = r • f /U  , 

Q ^ 


(13.14) 


where  we  have  employed  the  definition 


B = 2X  - r • B 

O 


(13.15) 


It  should  be  noted  tliat  Eqs.  (13.9)  and  (13.14)  are  exactly  the  saji\e  as  in 
classical  linear  elasticity. 

We  must  now  eliminate  H and  cp  from  the  representation  in  order  to  express 
Eq.  (13.5)  in  terms  of  B,  B , f and  V.  To  this  end  we  first  take  the  curl  of 
(13.3)  and  employ  (13. 3)^  ^ to  obtain 

V^H  = - V X R - a,V^V.  (13.16) 

We  now  siibstitute  from  (13.15)  into  (13.11)  to  obtain 


cp  = [ (r  • B + B )/2a,  ] - (otjoi  ) \ji  . 

’ ev  rv  O J.  Z J. 


(13.17) 


Taking  the  divergence  and  the  curl^  respectively,  of  (lJ.5)  and  employing 


1 3^  "ha  170 

' ^2,4^  


V [ (P,  +P,,)V''cp  + (b  + 2b„)y  ■ii  - a, ill]  + V • f ' ' = 0 , 
v^[|  + (b„  - b,)v^y  - a,y  ! - y x f = o . 


(13.18) 


(13.19) 


'2  "3 ' ' 4, 

Now,  substituting  from  (13.16)  and  (13.17)  into  (13.18)  and  (13.19)  and  re- 
arranging terms,  we  find 

■}  'I  A *\  /^  \ 

(13.20) 


(1-  = H (r  • B + B ) + V • f . 

X X'^~0  X' 


(1  - i^V^)V=  - K.,v^vx  B - vx?  Va,  , 


(13.21) 


where 


^ 2b^ 


a^^(X.  + 2iJ.)  ’ 


" ToTlhTa^'  ' 


j2  - **2 '^3 
*1 


2 

"2  - 2.^ 


(13.22) 


The  substitution  of  (13.3)  and  the  gradient  of  (13.17)  into  (13.3)^  yields 


U = B-Q''V(r’B+B  ) r— — ^ Vf  - — V X V , 

.\J  l'^  O A ■(■  2(il  2|^‘  ry^  ^ 


(13.23) 


where 


= (A  + ui,)/2(A  + 2m- ) 


(13.24) 


Thus,  finally  the  representation  consists  of  the  differential  equations  (13.9), 
(13.14),  (13.20)  and  (13,21)  along  with  the  expressions  (13.23)  and  (13,3)2. 

The  functions  B and  B reduce  to  Pap)covitch  functions  when  p.  and  p.  vanish. 

O 12 


14 , Concentrated  Forces 

In  this  section  we  consider  first  the  conc.’.ntrated  force  and  then  the 

concentrated  relative  force  located  at  the  origin  in  an  infinite  isotropic  two- 

constitueiit  composite  medium.  In  an  infinite  medium  acted  upon  statically  by 

'■  (1) 

body  forces  f and  relative  body  forces  f we  have  Eqs.  (13.9),  (13.14),  (13.20). 
(13.3)  and  (13.21)  for  the  potential  functions  B,  B , i|i  and  V,  which  enable 
the  determination  of  the  displacement  fields  u and  w through  (13.23)  and  (13,3)_. 

In  the  case  of  the  concentrated  force  we  have  everyvjhere  and  ^ = 0 

outside  a region  v'  encompassin  ' the  origin  and  containing  a nonvanishing  field 
of  parallel  forces  f.  A concentrated  force  is  defined  in  the  usual  way  by 


= lim  f 
v'-o  « 


f dV  . 


(14.1) 


52. 


For  the  case  uiider  consideration  Eqs. (13.9)^  (13.14)  and  (13.3)^  take  the  form 

shown  and  Eqs.  (13.20)  and  (13.21)  may  be  integrated  twice  to  give 


<1  - £^V^)\|(  = K.  V^(r  • B+B  ) , 

1 j.  /NW  ^ O ' 


(14.2) 


(i  - Jt^V^)V  = - 

2.  ^ £r~‘  ^ 


(14.3) 


Since,  for  the  infinite  medium,  solutions  of  equations  having  the  respective 


forms 


2 2 2 ~ 

V e = a , (1  - X V )9  = a , 


(14.4) 


can  be  written  in  the  respective  forms 


e - - i j dV  , "8  = — ^ I a (Q)  dV  , 


(14.5) 


r^  = [(X  - §)^  + (y  - Tj)^  + (z  - , 


(14.6) 


is  tlie  distance  from  the  field  point  p at  R to  the  source  point  Q at  r'  and 


2 2 2 2 2 2 
P • R - X + y + z , r ■ • r ' = I + Tl^  + r , 


(14.7) 


we  have  from  (13.9),  (13.14),  (13.2.0)  and  (13.21) 


. f dV 


^ ~ _ _1^  ’ ~ . 

4ny:  Jr,  •’  ®o  ~ 4m  J r 


r'  • f dV 


V 1 


V 1 


(14.8) 


^1 

4ttX^  V ^1  ® 


(r ' • B + B ) dV  , 

O ~ r\  ^ 


’^2 

V T V-  X B dV  . 

~ .2  r,  ~Q  ~ 

4nX^  V 1 


(14.9) 


(14.10) 


Then,  from  (14.1)  and  (14. 8)  in  the  usual  way,  we  obtain 


B = P/4mR , B =0 


(14.11) 


53. 


since 


lim  r = R,  lim  = 0 
V'-'O  ^ V'-*0~ 


(14.12) 


The  relations  in  (14.11)  naturally  are  the  same  as  in  the  classical  theory  of 
isotropic  linear  elasticity.  Since 


P = Pe  , 


(14,13) 


where  e denotes  the  unit  base  vector  in  the  z-direction,  we  have 
~z  ^ 


/ 2 2 2 
R • B = ZP/4T1VJ-  7r  +2  , V X B = (P/4TTijiR  )e„  , 


(14.14) 


where  e„  denotes  the  unit  base  vector  in  the  cylindrical  coordinate  0-direction 

A/D 

r is  the  magnitude  of  the  cylindrical  coordinate  radial  position  and,  of 


courESj  hsvs 


2 2 

r + z = R . 


(14.15) 


Substituting  from  (14,14)^  into  (14.9)  and  converting  to  cylindrical  coordinates, 
we  have 


„ « 2rr 

HP  nr" 

♦'-rr-;?  IJ  J ~ . ,z, j 

8n  oo  ) 


(14.16) 


where  in  cylindrical  coordinates 


r^  = [r^  + r'^-2r?  cos(0-0')+  (z  ’ 


(14.17) 


and 


-2  2 2 
r = X 


+ Y^,  = 1^+11^,  0 =tan  ^ x/y  , 0'  =tan~^  5/"n  . (14.18) 


At  this  point  it  should  be  noted  that  since  the  integral  over  6''  in  (14.16)  has 
an  interval  of  2rr  which  is  the  period  of  cos(0-G')  in  (14.17),  the  resulting 
expression  for  'll  in  (14.16)  is  actually  independent  of  0. 
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Since  V X B is  spherically  symmetric  it  is  advantageous  tc  return  to 
(14.13)  rather  than  to  use  (14.10)  directly  to  obtain  V because  cm  ordinary 
differential  equation  in  R results.  Substituting  from  (14. 14)^  into  (14,3), 
we  obtain 

(1- 0^)V=  (K_P/4np,R^)e-  . (14.19) 

z ~ ^ u 


As  the  solution  of  (14,19)  we  can  take 


X = ''e  . 


<14.20) 


because 


V • V ==  0,  Sfep,/dR  » 0 


(14.21) 


Then  V takes  the  spherically  symmetric  form 


V“ = R~"d (R~d/dR)/dR  , 


(14.22) 


and  substituting  from  (14.20)  emd  (14.22)  into  (14.19)  and  employing  (14,21)2, 


we  obtain 


[1  - X2R~^d(R^d/dR)/dR]  Vq  = H2P/4TpR^  . 


(14.23) 


On  account  of  the  relation 


-12  2-22 
^ j**  /f»TT  ^ _ r»  /^r>\  /ATi 

X\  \X\V  Q / / \i'.  VA  V0/  v*x\,y  y v**s,  y 


Eq. (14.23)  can  be  written 


RVg  - X^d^  (RVg)/dR^  =4  H2P/4T|iR, 


(14.25) 


the  inhomogeneous  solution  of  which  is 

R 


'^2*’  r 1 . , (R-  s)  , 

~ sinh  — 5 ds  . 

G 4np,X.,R  J s X- 

2 O 2 


(14.26) 


Tlius,  i‘  and  V for  the  concentrated  body  force  have  been  written  as  definite 
integrals,  and  we  carry  the  solution  for  the  concentrated  force  no  further. 
Clearly,  if  the  coupling  coefficients  and  P,,  vanish,  the  solution  reduces 
to  that  of  the  classical  theory  of  elasticity. 
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For  the  case  of  small  and  which  should  be  most  common,  asymptotic 
representations  of  the  solution  for  if  and  Vg  in  terms  of  simple  functions  for 
R « and  R » can  readily  be  obtained.  Ihese  asymptotic  solu- 

tions cam  be  matched  in  the  intermediate  region,  however,  although  Vg  can  be 
matched  relatively  easily  since  it  satisfies  an  ordinary  differential  equation, 
the  matching  of  ij.  requires  some  effort  because  it  satisfies  a partial  differ- 
ential equation.  Consequently,  these  asymptotic  representations  will  not  be 
treated  here. 

A 

In  the  case  of  the  concentrated  relative  force  Me  have  f = 0 everywhere 

A / T \ , , 

and  f = 0 outside  a region  V' enclosing  the  origin  and  containing  a non- 

A /I  \ 

vanishing  field  of  parallel  relative  forces  f ' . A concentrated  relative  force 

is  defined  by 

lim  f f^^^  dV.  (14.27) 

~ V'-O  ~ 

For  this  case  Eqs. (13.9)  and  (13.14)  take  the  form  shown  and  by  virtue  of 
(14,8),  we  have 

B ==  0 , B = 0 . (14.28) 

^ ^ o 

Tn  froin  (13.21).  we  have 


v^d  - = V • f , 

(14.29) 

v^d  - X^V^)V  = - Vxf^^Va,  . 

(14.30) 

since  for  the  infinite  medium  solutions  of  the  equation  having 

the  form 

2 2 2 
V (1  - X V^)0  = a , 

(14,31) 

48 

can  be  written  in  the  form 

1 -1 
^ J ^1  ^ ^ 

(14.32) 

V 
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we  have,  from  (14.29)  and  (14.30), 


(1)  -1  ”^1'^'^! 
* J f< 


(14.33) 


1 r -1 

V = n X [r  (1  - 

4na  ,1  ~ 1 

1 c 


e ‘^)f'  ']dS 


*4^  I 

J*  XT 


^ (1)  -1 

f'^^Xv[r.  (1-e  ))dV, 

~ ~ 1 


(14,34) 


and  we  note  that  the  surface  integrals  in  (14,33)  and  (14.34)  VcUiish  because 
= 0 outside  V'.  By  virtue  of  (14,12)  and  (14.27),  in  the  limit  V'  0 


Eqs.  (14.33)  and  (14.34),  respectively,  reduce  to 

■B/jJ 


1 rl  II  1 rl  2 i 

= (1-^  >]>  :x  = ^gxs[R  )]- 


(14.35) 


which,  with  (14.28),  are  the  stress  functions  for  the  concentrated  relative 
force . 
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